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HYDRODYNAMICS OF FLOW INTO CURB-OPENING INLETS © af in 


By Ri Richard J. . Wasley, 1A.M.ASCE 


4 


Gat 


curb- -opening pety is presented. ‘The investigation is divided into two coe 
(1) the instantaneous failure of a dam retaining a triangular reservoir super- 
Posed o on (2) a 


a certain incoming velocity distribution. ae verifica- 
is given 


at present, sev 


eral methods are available for c from road 
surfaces into storm drainage systems. The scope of this inital is confined to 
opening inlets without depression in the gutter line. 2 
The analysis concerns an open channel (Fig. on a ‘continuous 
ysis pen Pp 
grade; it is of triangular cross- section, one side of which is nearly vertical 
(the tw two o sides of the section are mutually perpendicular). At a distance from 
_ the upper er end of the channel sufficient to establish steady, uniform flow condi- a 
iy (see * ‘Comparison of Theory and Experiment”), the vertical side ends 
abruptly. “The termination of this wall marks the beginning of the curb inlet 
Note. —Discussion open until January 1, 1962, Separate discussions should be sub- 
mitted for the individual papers inthis symposium, To extend the closing date one month, 
_ a written request must be filed with the Executive Secretary, ASCE, This paper is part 7 7 


_ of the copyrighted Journal of the Engineering Mechanics Division, Proceedings spinel the | 
American Society of Civil Engineers, Vol. 87, No. EM 4, August, — 


1 Research Asst., Dept. of Civ, Engrg., Stanford Univ., Stanford, Calif, 
2 “Hydrodynamics of Flow into Curb-Opening Inlets,” 


by R, J, Wasley, Ph, D, Dis- — 

“Dept of Stanford Univ., Stanford, Calif., Nov., 
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1960 (available as Tech, Report No, 6, 
ept. of Civ, Engre., Stanford Univ., Stanford, Calif. 
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Section under study. The length of opening is s1 sufficient all the 
pass over the edge of the inlet lip. Rainfall and lateral inflow are not consid-— 

ta Slopes greater than the critical a are considered; this is easier to analyze 
‘theoretically than the subcritical case. In the latter situation, the complication 
7 toward the opening, since the transverse . slope is almost always ‘supercritical. d 
7 In ieee a certain amount of drawdown will occur before the control section 


intothe section, makes the in direction, and passes over 
7 the edge of the inlet lip. The one-dimensional component velocities ee 
by averaging the point velocities in the vertical direction) and the. depths, as 
: functions of position, are the parameters of interest and are documented theo- — 
the results are verified experimentally. 


THEORETICAL DE’ DEVELOPMENT 
Notation. _—The symbols adopted for use i in this paper are defined 
where they first appear and are arranged alphabetically, for convenience of 
_ General Problem.—The water passes down the channel ina uniform state 
. as given | by the flow rate, Q, and the pertinent one-dimensional velocity com- 
_ ponent, u. (Fig. 1) shows the essential nature of the problem and the relevant — 
parameters. Upon reaching the section ODA the water begins the transition in — 
direction in order to pass over the inlet lip, OC. The general movement of the 
_ fluid down the channel suggests considering the total flow as a main flow in the 
x-direction combined with a cross flow in the y- -direction. Such an approach | 
. resolves the actual flow into two simpler components, each of which may be 
treated theoretically. Hence the investigation is divided into two phases: (1) 7 
the main flow, i.e.,a steady open- channel flow witha certain incoming ny 
distribution (either experimental or theoretical), upon which is superposed (2) 
the cross flow, which corresponds approximately to the instantaneous failure > 
of a dam, behind which there is a triangular reservoir. 
a By this procedure (i.e., by initially considering the inlet problem as actually 
two problems, solving them separately, and later combining the results),a dif- 
. ficult steady flow situation is transformed into a relatively simpler problem 
- involving an analysis of an unsteady flow coupled with a steady flow. There is 
a similarity in approach here to the slender-body method utilized frequently i. 
in aerodynamics in which two vector fields are superposed. present prob-| 


lem, however, is more difficult because of the rotationality of the main flowin | 


the transverse direction, 
ex. From Fig. 2, it may be 5 seen that because ‘of the assumption of supercritical 


_ flow in the downstream direction, there is no Vv component at time t=0 (x=0), 
: _ and hence the dam failure is directly applicable without considering some ini- 

a " Some interaction between the two phases involved in n this superposition is 


unaccounted for, but not so much as might be originally expected (as seen by 


— 
— 
| 
a 
— 
— 
4 
— 
4 
| 
4 
— 


2 TIVNOLLOGS CNV NVId MOTA OMLVINSHOS—" 


900 NOILO3SS 


NN 3YNSNI O 


fa 
Bn 
og 


> 


\ 
: 
= 


INLETS 

q 
A 

— 
— 2 


August, 1961 

coefficients except the friction om experimental information fog 

handbooks) an average friction | factor ‘across the a can be given | by3 


PERI 

=4R is th the. e “equivalent Ri is hydraulic ra radius, 
and u is the average velocity of the channel; these se quantities c can be readily 

_ Considering a certain point on this transverse section, ar and rewriting scttadel 


e uation, one finds that. 

in which Co with the depth being the 


radius. 

sucha s simple ‘representation of the distribution as Eq. 2, it could 
_ be expected that complete agreement with experiment would be more coinci- 
dence than consequence of the theory. It is seen, however, from Fig. 12 12 that “4 
in general the distribution is quite well defined. 


It is noted that Eq. 2 does not account for the curb shear which would cause 
_ the velocity there to approach zero. ‘If the equation could be modified to ac- 
count for this, only a small correction would result. a 
Some arguments could be raised that an “average” friction factor, defined 
by Eq. 1, might be questioned because of the complicated variation of f f across” 
channel. Also, the use of the “equivalent diameter” (d= 4R) for such 


shallow channel might cause some concern, 
Instantaneous Failure of Dam. —It is readily possible to derive the differ- 
ential equations governing the motion of fluid down a channel after instantane-_ 


failure a dam. J. J. ‘Stoker4 has two- ‘dimensional linearized 


will be focused upon a one- -dimensional approach. 
: _ Considering this phase of the total problem alone, a sketch (Fig. 3) can be 
| 7 drawn indicating the parameters necessary to describe the dam break problem. 
°1 (A variable with an asterisk superscript refers to — measured to the 
For ease of experimental measurement, the y- and are de- 
: fined as in Fig. 2 and shown dotted in Fig. 3. However, owing to the small 
pe 9, the axes may be rotated slightly counter-clockwise as shown, to aid ~ 


in the mathematical treatment. The slope, of course, is still considered in the 
analytical solution. _ The error introduced when comparing the solution with — 
_ experimental results using the nonrotated coordinate system is negligible. 
(Actually however, the vertical depths, given by the distance, d (see Fig. 1), 
were the measurements taken. ) Note that the channel bottom is assumed to 

continue in the negative y- -direction. 


38 “Elementary Fluid Mechanics, by. J. K, Vennard, John John Wiley and Sons, , Inc., 


“Water Wave Waves,” Stoker, Interscience Publishers, Inc., New York, N. Y., 
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‘angle approximations, i. , sin 80 and sin ~ 


Here | is is the tection slope (defined this” term is included in the 


of the governing equations but in the solution. By 


those terms i in w differential of second onder 


q 


vig 


Mi 


which is the > appropriate ite form of the dynamic equation for the problem consi- 


10 


on 


which is the desired form of the continuity equation for the present problem. 
ns, The exact solution to Eq. 3 (excluding 6¢) and Eq. 4 is rather = 
but by a technique proposed by R. F. Dressler,° this can be achieved. Dres- 
sler states that in general the exact solution to every non-degenerate unsteady 
-water-wave problem (governed by the above equations) in a straight channel — 
at an arbitrary slope canbe obtained in terms of the complete elliptic ean 


- of the second kind. By means of a non-Newtonian reference frame (i.e., an 


— 


“Unsteady non-linear waves in sloping channels,” by R. F, Dressler, Proceedings, 


Royal Soc, of London, Vol, 229A, March 1958, pp. 186- -198, ae iat 
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"accelerating coordinate e system) every y such wave problem can be : replaced by 
an associated problem for a horizontal channel. Then the equations are trans- 
_ formed into the Euler-Poisson equation through the hodograph inversion meth- : 
od. _ Application of the Riemann method of integration yields an auxiliary a 
— tion which is transformed into a Legendre function and then into the — 7 
- integral. The resulting plots given by Dressler enable the depth and — 
a mensional velocity profiles at various times to be constructed. ee . 


os The parameters can be non- -dimensionalized by the following crea 


“The oneevae are self- -explanatory | in the above cases. Several illustrative a 


a ie (in non- “dimensional form) are given in Fig. 5 for times of 0, 0. 3, 2. 6, - 


| 


FIG, 4, ._—DERIVATION SKETCH 


Superposition. —There isa natural breakdown of the total curb-opening inlet 
flow into two definite sections, labelled Zone I and Zone II, as indicated in Fig. 
6. In the flow field a definite curve (OB) is formed at which the water appears — - 
to be “sloughing off” and forming a crest. This curve is the boundary in the y y- 
e_ of the negative surge crest as it passes downstream due to the x- © 
component of velocity. In other words, the information that the curb has ter- 
-minated at the origin 0, is transmitted in the y y-direction at a certain finite 7 
which shall be derived. Superposed on is a motion in the downstream 
‘direction at a rate governed by the incoming velocity distribution. The delin- 
= of OB was visually observable inthe experimental phase of the program > 7 
a = for certain configurations (4 equals 1 on 24 and 1 on 12 for all 0) » al-— 
though it was measurable by instrument for all positions investigated. 
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<a Before Zones I or II can be determined, it is necessary ‘to Cotermine the i. 
movement of the crest of the negative surge inthe transverse direction. ‘Know- 


= the velocity of propagation of a small ee c= = 


im 


where in this « case, the asterisk superscripts refer to ) measurements taken | to 
the boundary « of tone I (Fig. 6). Substitution of Eq. 7 into the above leads to 


ion between the e limits of 0 and and. and results. 


7 and 9 give re the param — the line OB. the 
Be tis ) etisninated from these two equations a linear relationship in “cal 


4 
integrating between the limits of 0 and 
where t < . By non-dimensionalizing, according to Eqs. 5,there results 
<2, 
Zone I.—Now re = imensional absolute 
locity, U, in terms of the equivalent differential quantity, one finds that 
im 
4 
. ® iim” 
— 


. The first superposition of the analy 


‘In Zone I no change in upstream conditions is possible (for mpeverttien 
conditions), and the surface level and velocities at various points are easily 
be found for continuing the solution into the interior of Zone II. . One way = 
a this involves the description of the trajectory of a given fluid element. 
As it moves transversely across the channel it simultaneously translates in 


the longitudinal direction according to the of the velocity dis- 


“FIG. —TRANSLATED COORDINATE SYSTEM 


tribution. Combination of these two motions determines the complete velocity, 
_ The analytical definition of the limiting streamline in the zone under study 
is one of the more important parts of the solution. By the relatively simple 
result derived, the clear length of opening necessary for complete interception _ 
can be obtained immediately, knowing the upstream parameters of the flow. 
i. From Dressler’ 8 exact solution of the hydraulic problem é associated with 
the instantaneous dam failure, the last particle of flow is 3 governed by a very 
simple equation influenced by gravity alone, that is, frictional and surface ten-— 


effects are neglected. At t = ), 0 one has and x = Co 


y= Yo» 
8, and the motion of this particle is given — 7 ee ee 


=-g9 


linear translation of the original system. 


— 
on Hence, it is a very simple mat — 
4 downstream boundary of Zone I 
inition for non-dimensionalizing the x-coordinate depends upon the theoretical _ 
AE 
= 
x oO 
— 
jj. To simplify the solution of Eq sill 
— system as shown in Fig.7froma [im 
= 
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The solution is readily obtained 


as (0 y 


om 
Substituting 13 into Eq. 14, and integrating the limits of 0a 


er 2, 5, the 


which isa linear in x’ y’. terms of the original 


Hence Eq. 12 can be rewritten as, 
using the condition that when t' = 0, y' =y' =0. 
4 Eq. 2, using the modified reference frame, with the relationship z = 
it is found that — 
= 
; 
‘7 
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~ Xotheo. 
TABLE .—BOUNE NDARY EQUATIONS FOR ZONES I AND 


: y = A, 0 < x< 


< 


Vo 


Xotheo. = 


_ Acomment may be made here which applies to all the streamlines in Zone 

Bi nine the boundary, since friction is not considered. As depth decreases, 

the velocity decreases in uniform, steady flow. It would seem, therefore, that 
some adjustment should be made for the downstream component of velocity. © 
However, it is unlikely that the fluid would decrease immediately to uniform 
velocity (given by Manning’ s or Darcy’s equation) for each incremental 

distance inthe region under study; inertial effects would have to be considered. 
Hence in the short Comanees involved, the particles are assumed not to de- ; 

crease significantly | in velocity in the x-direction due to decrease in depth. 7 
_ The delineation of the various intermediate streamlines in the interior of. 
Zone II is unfortunately not so straightforward as the analysis presented for 
the boundaries of the two zones and the interior lal Zone I, and a finite incre- 
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initial position 
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n Fig. 5 at 


Consider a t a transverse at an arbitrary time, ¢ (Fig. 3). ‘Assume the 
de pths and velocities (in non-dimensional form) at ve various points as having 
ane previously been obtained. A particle, Pj, located on OB, is then selected. This” 

corresponds to any given desired streamline in the x-direction in Zone I de- 
pending upon the assumed, example, the streamline dividing the flow 


in halt 


‘Now. a certain interval of time, 68, is assumed to have elapsed. The 2 veloc- 
ity Vy, at time ¢ + 
5. _ This allows an average of the two seheaibing Vip and V1) to be taken, thus 
_ determining approximately the position of the particle. This process is con- _ 
; tinued until the position of Pj has stabilized for the assumed small interval of _ 
time. ‘Two cycles are sufficient for practical application of the method. — 
The actual final position (in terms of & and 9) that the particular particl 
takes: ‘in 54 can be found quite readily by determining | the average down- 
stream 1m velocity, Uavg. » the time interval. 


— 
Rip | Sip | Zip (8) 
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ime for 


ad time, 


Fi 


Irom 


k (k = 2/(2 + V2) appears; this is because of the nature of the non- ~dimensional - 
desired, it is possible to describe adequately any general experimental 
incoming velocity distribution, including the fact that the velocity approaches — 
zero at the curb face, by means of suitable power expressions with arbitrary 7 


j exponents. | ‘By a ‘similar procedure to that developed previously, power series a 


The iteration procedure is given in Table 2. “Under the column 6X, a factor 


expressions can be obtained which determine the boundaries of Zones I and II. -_ 


_ (In this case, closed form solutions cannot be obtained, and the pertinent dif-_ 
ferential equations are expanded by use of the binomial theorem and then inte- 


7 grated term ads term in the convergent domain). 6 
COMPARISON OF THEORY AND EXPERIMENT 
a A matrix of apparatus configurations was investigated experimentally. This 7 
information is most easily and clearly presented in a plot such as Fig. 8. The 
solid circles indicate that the measurements were taken for these arrange- 
a It is seen from this figure that the lower limit of supercritical flow possibly 
occurred between the longitudinal slopes? of 1/2% and 1%. Since the theory | 
was initially predicated upon uniform and super- -critical flow - being maintained 


_ A dimensionless plan of the inlet section is presented in Fig. 9. Zones I and 
I are delineated from the analytical methods previously developed, and from 


i theory certain intermediate streamlines (at 25%, 50%, 75% and 90% of the tot: tal 


“Hydrodynamics of Flow ‘into Curb-Opening Inlets, by R, Wasley, Ph, Dis- 
sertation, Stanford Univ., Stanford, Calif., No., 1960. 
7 “Hydrodynamics of Flow into Curb-Opening Inlets,” by R, J. Wasley, Ph, D, Dis- 


-sertation, Stanford U Univ., ‘Stanford, Calif., Nov., 1960, 


(43) @4) | (as) a6) | a7) | 8) | | (20), (21) (22%) | (23) 
a 
a 
iia 
— 
g ig 
| 
— 


flow) are plotted. The experimental data d defining the boundaries o of the flow 
are shown, with the exception of the data for the 1 on 96 transverse slopes; : 

this slope introduced additional phenomena which invalidated the adequate rep- 
resentation of conditions by the present theory. ar shallow depths with small 


friction and surface tension to become r more 

q ee The agreement between experiment and theory is quite good. i The experi- | 
— mental intermediate streamlines, especially the 50% streamline, corresponding 7 


the divisions above, ‘compare with theory. The depths 


© 300% =! 24, 


FIG, 9.—STUDY SE SECTION PLAN—THEORETICAL AND EXPERIMENTAL RESULTS 


and velocities on the various experimental streamlines do not coincide with 
theory quite a: as well, but they a are still fairly satisfactory, particularly in in —_ 


formulas necessary to convert the non-dimensional theoretical cae in Fig. 7 
9 to the actual corresponding physical values are — in the plot ‘itself a1 and 
is seen that the theoretical boundary a sudden change 
in direction at t point B. However, it is clear fromthe fact that the 90% stream- 
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FIG, 10,—DISTRIBUTION OF FLOW ALONG LIP—THEORETICAL 
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line is relatively a considerable ee eta this boundary, that no practical 
significance should be attached to this. 
_ It is desirable for practical purposes to determine the distribution of flow 
at the lip of the inlet in order to locate the position where some certain per- 
centage of flow has been intercepted. This" distribution is ‘plotted in Fig. 
(for ¥ = 0) and ‘compared with the experimental results as shown (Q = 0. 725° 
efs, = 3.00%, and = 1:24). The variation is almost linear - both theoret- 
ically and experimentally. a situation could conceivably be used advan- 
tageously. If merely a good estimate of this quantity is desired, then to con-— 
"sider this an exact linear velationship (knowing Xo) wo would be more than ade- 
, 2 was stated that one of the assumptions 1 used in the development « of the an- 
_— was the consideration of uniform flow upstream from the re- 


50 


ane for this particular point were essentially identical. For all eositions, 

_ configurations, and flow rates investigated, this situation existed, thus indica-— 
ing that the uniform flow requirement was satisfied. 
It was mentioned that Eq. 2, although quite simple, defined the actunt exper- 
imental velocity distribution adequately. A typical transverse distribution is - 
given in Fig. 12 for Q = 0. 25 cfs, $0 = = 3.00%, and % = 1:12, which — 
surface fluctuation of high frequency (probably caused by turbulence effects) 
was noticed in all runs made; addition , surface tension 


— 


: 
- tions (x = 0 ft, -1.00 ft, -3.00 ft, -5.00 ft, and -7.00 ft) for a given configuration - ae 
| 
ae (% = 1.00%, 6 = 1:12, and y = 0.50 ft) and flow rate (Q = 2.98 cfs) is shown — ae 
3 
‘ 
: 
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ine were observable, particularly in the shallower transverse slope ap- 
paratus configurations. However, no problems were encountered in the actual 
maeaourement of depths or velocities « of flow. It was relatively easy | to repeat 4 
‘readings to the limit of the accuracy (0.001 ft) of the depth-measuring point 
gage. For velocity determination, Pitot probe results were duplicable to a J 
in about 2% in most cases | (the most critical situation having a maximum error 
work covered a sufficient. range of slopes (both transverse and longi- 
tudinal) to describe adequately most practical situations. It is quite certain 
that the theory developed in this paper is applicable to longitudinal slopes 
greater than 5%. It is also evident that the analysis is satisfactory for $o = : 
0. .5% (see Fig. 9), but below this slope, the theoretical results may be inac- " 4 
7 urate. The steeper cross-slopes were best when compared with theory; this 
reasonable because of the assumptions in the analytical approach. 
_ The range of flow rates investigated was from 0.0032 to 2. 98 cfs. This cov- 
ers the practical variation of flows in most situations. The quantity of flow 
does not invalidate the theory except if it becomes so small as to decrease the 
velocity to the point where frictional and surface tension effects become ‘signif- 
a The ery is ; simple anc and straightforward, , depending only on knowledge of 
the upstream parameters ‘~ Go, Bo, and f in its use, and i is quite satisfactory 


tations) from Figs. and 10. 


direction this research was conducted. "Discussions with Carl F. Izzard, Chief, 
Division of Hydraulic Research, ‘Bureau of Public Roads, provided much useful — 
_backgroundinformation. The assistance of F. G. Baker, J. Duncan, R. L. Street, " 

and Y. Wang is also acknowledged. Financial | support for this research was” 
the S. Department of Commerce, of Public Roads. 5 
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a defined by yin per sec; 
: = absolute celerity of wave in still v water, in feet per second, de- 


= vertical depth of water, in tt, defined by cos cos 


‘the friction factor; 


ve 
as 
3 
a : q formation on the depths and velocities (with directions) in the inlet section and 
— ‘tial interception compu- wea 
| ¥ 
— | 
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an 
= rate of flow, in cfs; ae 


the hydr raulic iain, in ft; "a 

time, insec; 


the x- and y- “components of velocity | (average in vertical ), re- 


spectively, in fps; 
= the x-component of velocity (average o of section), in — 


the coordinate axes and distances measured ‘parallel thereto; 


max 
al pecific wei ht of wa water, in lb per cu ft; Pe s <i 
the fri the friction on slope, , defined by ss Bg ss a, 4 


the density of of water, in cu f 


fo the slope, in 
( to crest or line OB; 
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IN’ THE PRESENCE OF PRIMARY BENDING MOMENTS 


F. Masur, ! I, Chang,” and L. 
we ys 


A systematic > method of ‘euslyaing the stability « of frames against - sidesway 
= is developed. Two possible avenues of approach are presented: (1)a_ 
slope- deflection equation method based on equilibrium c considerations; and (2) ‘ 
a moment distribution procedure, These techniques are well-known in connec- 
‘tion with the analysis of the stability of frames subjected to loads at the joints» 
only; the present paper represents an extension to 1 the case ir in which 
primary bending moments occur prior to buckling. — 


if a framed § structure of the e type shown in n Fig. 1 is subjected to increasing : 


diagram. Collapse is said tooccur when that branchof the diagram 
has been reached on which the load magnitude decreases with increasing | de-— 
ee. linear, elastic, first- approximation analysis of the structural henge 


_ does” not exhibit such a maximum, Actual collapse is is due to factors not in- 


4 
_ Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
A written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 87, No. EM 4, August, 1961. _ 
ol Prof. of Engrg. Mechanics, Univ. of Michigan, Ann Arbor, , Mich. + ; 
2 Stress Analyst, Scientific Design Co., New York, N.Y. 9 ; 
3 Senior Research Scientist, Inst. of Science and Tech., Univ. of M Michigan, Ann Arbor, 
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_ cluded in this type of analysis; among these are chiefly inelasticity and the 
second order effects associated with “buckling” instability, or a combination — “4 


of both. These combined effects are difficult to analyze, although some pro- 
gress | has been made in recent years.” 4 The usual approach | is to. assess the 


‘The present paper is concerned with elastic instability alone. The problem | 


has been investigated extensively in the past; however, in almost all previous | a 
investigations of the subject, the actual loading system consists of, or is re- 
_ placed by one consisting of, forces that are applied at the joints —_. —e un- 


buckled structure : is therefore assumed to be in its virginal state, that is, its 
members are straight and no “primary” bending moments are present. The 
replacement of the actual loading system by one that is (for each member) 
_ tatically equivalent toit is usually justified by 1 the assumption that very | small 
errors are thus introduced. Pessie assumption is fortified by a result of E. 
"4 “The Stability of Tall Buildings,” by R. H. Wood, Proceedings, ICE, Vol. 11, 1958, 


_p. 69; discussion by M. R. Horne, et al., of “The Stability of Tall Buildings,” ws R. A. 
Wood, Vol. 12,1959, p. 502,00 
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FRAMES 


_— Chwalla5,6 6 who has made an exact study of an example, including | the effect of 7 
primary bending moments, Unfortunately, Chwalla’s work is almost the — 
one to date of which the writers are aware; moreover, it is unsystematic 2 and 
E Chilver,? A schematic discussion of the problem has been presented by AL 
WL Chilver,? who replaces the actual structure by a much simpler one thought _ 
to be equivalent. The present study seems to confirm the relatively insignifi- " 

Se effect of primary bending moments onthe elastic stability of structures of 7 

A the type shown in Fig. 1. _Nevertheless, the results presented herein are con- 

7 _ sidered to be of potential value, in the sense that they may represent a first 7 
"step toward a more comprehensive plastic instability analysis, for which pri- > 
mary” bending moments may not be considered negligible a priori. Similar 
doubts, as applied to the case, have recently been expressed by 

‘ aes Two avenues of approach are presented. These are: (1) an equilibrium 
analysis leading toa eet of “slope- deflection equations” and (2) a moment- 
distribution | technique. Both of these methods can be seen to be extensions of 

_ classical methods that do not take the effect of primary bending moments into © 

account and that are treated comprehensively elsewhere. 6 Of course these 
two formulations are not the only ones possible, For example, in n analogy with 

‘more c conventional problems of structural analysis, a solution may be obtained 

through the use of the conditions of compatible deformations. This type _ 

procedure is described more fully elsewhere, but, together with other 


ble lines of has been omitted here for sake ‘of i! 


‘SLOPE- METHOD 
- In het pte: a brief resume is given of the development of the standard 
slope-deflection equations. Although this subject is treated exhaustively in 
P ‘standard treatises, 6 a cursory summary is considered necessary in n order to 
make the present discussion reasonably self-contained, 
_ A typical member connecting joints i and j is shown in Fig. 2. It is sub- 
: jected to an 1 axial force Pjj, which is considered positive in compression. It 
is also subjected to a lateral load qij (x), which is orange if its direction is 
_ clockwise relative to the x-axis pointing from ito j. The rotation of joint i is 
If Pij and Bij are, respectively, the bar rotation and the slope of the 
‘flected beam at joint i relative to the deflected cord (i, j), then the continuity 
_ of the structure requires that 


‘if f all terms are defined to be positive if clockwise. 


“Die Stabilitat lotrecht belasteter Rechteckramen, ” by! E. Chwalla, Der Bauingen- 


ieur, 
Inc., 


6 Buckling Strength of Metal Structures, ed F. Bleich, McGraw-Hill Book Co., 


7 “Buckling of a Simple Portal Frame, ” by A. H. Chilver, Journal of the Mechanics — 
of Physical Solids, Vol. 5,1956,p.18. 
Discussion by Eugene Guillard of “Influence of Partial Base Fixity o on ‘Frame Sta- 

bility,” by T. V. Galambos, Proceedin s, ASCE, Vol. 86, No. ST 5, May, 1960, p. 153. *, 

ih. “On the Buckling Strength of Frames,” by I. C. Chang, thesis presented to the Illin- 

Inst. of Tech. -» in Chicago, Ill., in 1958, in partial fulfillment of the requirements 
for the degree of Philosophy. 
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| ening momest (also” ‘positive clockwise) pers. on the bar (i, j) a at 


4 


in have been retained their omission may to. 
ambiguities, If the bar in question is prismatic, that is, if its modulus of elas-— 


4 


In these equations K duneniehe on stiffness cima while C may be in- 
= as the “ carry-over factor.” Both are tabulated a as s functions of the 


Fie. 2. —TYPICAL MEMBER 


‘there | is s slightly different from the present notation, For the case of tension, 
> ‘that is, when P is negative, 9 becomes imaginary and the trigonometric func- 
tions in Eqs. 3 are replaced by their hyperbolic equivalents. For 9 =0 the 
"values of K and Cc approach the conventional 4 EI/L and 1/2, respectively. 


Ri 10 ’ “Extended Tables of Stiffness and Carry-Over Factor for Structural Members 
Under Axial Load,” by E. E. Lundquist and W.D. Kroll, Wartime Report L-255, N.A.C. A. 


$ — 
— 
— 
4 
; 

ahs, 
as 
| 

— 


In Eq. 2. the ter term rm Mf represents a a “fixed-end moment. ” This c: can be ex- 
pressed 


(4a) 


(a) =a Yo = = ‘constant: 


te 
L 
1 5 
aconcentrated force of Q at adistance o L from joint 
moment equations of equilibrium at joint i that 


in which the summation extends over all the bars connected at joint se The 
- eee of equations | so . obtained is equal to the number of unknown joint To- 
_ tations There is | an equation of equilibrium against ‘rotation for each 
number of equations so obtained equals the number of bar shears, “Finally, 


two force equations of © equilibrium per joint correspond to the number of un- 
mown axial bar forces plus the number of independent bar rotations that the 
‘Kinematic tonstraints of the structure permit. = 
; In the present paper this is applied to the example shown in Fig. 1; this ex- 
ample has been chosen in order to afford a comparison between the method of 
: analysis described in the present paper and the numerical results i 
elsewhere. 5 Complete symmetry | of both the structure and the loading system 
_ From the equations of equilibrium of the structure as a — it — 
‘that the axial fo1 forces in the bars are given nby 
Pap = 


In] Eas. ' 
‘positive t to the right. From the symmetry of the problem w = 0 in the present | 


jeer suede, the term is included in Eqs. 7 for the purpose of later 
comparison with nonsymmetric configurations, The redundant horizontal re- 


— 
| k q a 
a 
if 
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a The enn of the le left column requires that 


in which represent w=0: in the present case 
Eqs. 8 and 9 are combined to eliminate the joint rotation @ Be When Eqs. 3, _ 
4, 5(b), and 7 are substituted this leads to the following relationship between 


ia 

92 

= 


SS Eq. 10(d) H can be determined explicitly in| terms of Q. Br, ae 


i- 


Up to this point the analysis is conventional and agrees with 


cussions given elsewhere. 5,6 The question arises as to whether the equilib- 

4 rium configuration si so obtained is stable or not. In the present paper, instability 
is postulated if, in addition to the configuration just _ obtained, there also exists. 
one which is adjacent to it while the loads retain their magnitude and direction 

and their position relative to the moving frame. (This definition of 


Euler and has been discussed exhaustively by H. Ziegler,11 who has recon- 
firmed its validity as well as its limitations. r ——— of such a condi-- 


implies one of two possibilitie 


(1) the: load deflection curve has a with 


. my (2) there exists a bifurcation of the load- .deflection relationship, which, ir in 
present example, signifies the possibility of impending sidesway. 
_ The first case, although of potential interest in problems involving sym- 


metry as in the present example, is usually of primary ieupertones only if such 


> 


Vol. 4, “Academic Press, New York, 1956, P. 351, 


ew of Eqs. 2, t L <4 
— 
ee eee ee (10b ind 
: 
ia of the columns and of 
nthe Cross beam, respectively. Eds. De Simpinied greatly by consider- 
ing that is usually a very small quantity. When the appropriate trigonometric | 
g y 
, 4 expressions are expanded ina Taylor series near W = 0, and only the lowest a 
a terms are included, this leads to the approximate relation == mz 
4 
— 
Advances in Applied Mechan- 


symmetry is not the example under 
= problem of general interest is that of the bifurcation of the equilibrium § 
- configuration, when an asymmetric configuration (sidesway) is imposed on 
the existing state of the structure; this occurs when the load Q reaches its | 
For this adjacent state » of equilibrium Eqs. 2 and 6 continue i be valid; of 
_ however, the moments M are replaced by M + AM, the rotations p by ptAp, 
_ and so forth, Furthermore , the ¢ coefficients K and C also change and must be q 
7 replaced by K+ AK and C + AC, as 1s will be demonstrated later. if these new = 
aA terms are then inserted in Eqs. 2 and 6, if the original Eqs. 2 and 6 are sub- _ 
tracted from the ensuing relationships, and if all terms of order higher than _— 
the first in the incremental terms are deleted, the incremental equations of — 
moment equilibrium take the form == | an. 


and 
4 Eqs. 12, together with appropriate incremental equations of bar equilibrium, — 
represent the > condition of neutral equilibrium | as outlined above. The terms > 
AG account for the effect of the prebuckling deformations and, hence, of the 
primary bending moments; they vanish if these effects are ignored, in which 


case the analysis becomes ‘identical with the analysis proposed in the standard — 


In Eqs. 12 the terms AK, AM, for "example, represent the change in the 
stiffness factor, fixed-end bending moment, and so f forth, associated with the a 
change in the axial force in the individual bars resulting from the < —_ a 
equilibrium configuration. That is, after linearization, 
AK = K'AP...... 
| 


ats 
ll 


in which K' = dK/dP = = (dK/d¢ ) (do /dP); the expression for AF is maeeeied 
_ Expressions for K', Cc’, for example, appearing in Eqs. 13, are obtained 7" = 
_ differentiating Eqs. 3 with respect to Qo . When ¢ is eliminated from the re-— 


a sulting formulas through the use of 3, then, some 


2 KL 

0) 
K 


= 


= 

q 
— 4 
— a 
— 
— 
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— a 


‘Eqs. ‘ 
pears ii in Eqs. 12, can be obtained ‘similarly; this is done by “considering Bes. 
4 4 and by differentiating the expression for F with respect to 9. The result 7. 


pends of course on the individual loading case, but the effect of this term is <7 : 
general small. It often vanishes altogether; this is, for example, true in the | 
Rettig case becai because, as will be shown later, A 4 P = 0 for the cross bar during © 


In most cases, the axial in the column is small compared with the 


Euler force associated witha simply supported column. With @ thus approach- 
ings zero, the enpreasions in Eqs. 14 14 may be appraximated, near 9 = 0, by 


+ 


a one end of the bar, for example, joint j, is hinged the equivalent it expres- 
“sions for Eqs. 2 and 12 become somewhat simpler and are given by a a 


which K’ given in 14 or 15. 


In view of Eqs. the terms” in Eqs. 12 AP repre-— 
senting the change in the axial force in the corresponding bar. This axial force ~ 
P may be found from the reaction if the equilibrium configuration of the struc- ‘-_ 


fore, functions of | the bar ‘Totations p as well as of the remaining (that | is, 


12 “Post-Buckling Strength of Redundant Trusses,” by E. F. Masur, 
ASCE, Vol. 119, 1954, p. 699. 


K 
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4 
— 

of 
: 
of the bar rotations p. In turn, some of these reactions may be established by = a 
of the annuationge of annilihrium for the whole ctructire: they sre there- 
| 


are used as independent variables the bar may be 


pressed inthe form © nd 
in which the summations over ‘the unknown rotations and redundant re- - 
i; In the present case, the kinematics of the structure shown in Fig. 1 —— 
that PAB= PCD = w/h, whereas PBC = °, The only redundant reaction in eon 


AB 


- in which the (ieath alialalis refers to the left column and the second equation to 
ts The incremental joint equations of equilibrium corresponding to Eq. ll are 


found similarly. In view of Eqs. 12, 13, 16, and 18 this leads eventually to the 


“AB Aw 


ay 


} 


A x, AH, and Aw. "he usual, a nontrivial solutionis possible if and only 

if the determinant of the coefficients vanishes; this the character- 
‘atte equation of neutral equilibrium, 
n the present example, the determinant splits into two 2 x 2) determin- 
virtue of the symmetry of the problem (Kap=K Kep 600 and 


| 
a _ Eq. 9 and from an equivalent equation “Pp. With W = ap 
this leads to 
(Ben * 48c * - 
...(20) 
which Eq. 20(a) refers to joint B and Eq. 20(b) to joint 
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0 forth). The two cases represented by these two subdeterminants are as 
(1) A 8 B= - 49¢ and Aw=0. This means that the symmetry of the con- 


| iguration is retained and the load- deflection diagram has reached its highest 7 
point. In most structural problems of practical significance the critical value 
Q so obtained is much higher than the one found from case; 


Fan), 
vanishing of the determinant of this 


Kap 


‘at 


(Kap + ac) 


in which r may be as a “moment distribution factor” at joint B as- 
- sociated with an antisymmetric deformation. The term @ B in Eqs. 22 is found 7 
from Eq. 8 and represents the effect of the primary bending moment for the 
example under consideration; if it is omitted, the analysis agrees econieatneatl de- 
P- If in Eqs. 22 the appropriate expressions s for t the » stiffness factors and their “a 


derivatives are substituted, and if 8 is the 


cot (1 - cot cot 6 
Eq must jointly w: with 10 governing ‘the value H. is 
interesting to note that, except for aslight difference in the second term, which 
_ accounts for the presence of a small axial force in the cross beam, Eq. 22(c) 
_ is identical with that given by Chwalla. 5 However, Chwalla’s work involves the > 
integration of several differential equations and the satisfaction of the asso- 
Clated boundary conditions, last term represents” the effect | of the pri- 
mary bending 1 moments because H = 0 if the e loads are applied at the joints, that 7 
Numerical values have been obtained fox, the special case I = = and L=h 
If a = 1/3 the critical force Q is given by 92 : = 1,775 as against o2 = 1,816 for 
a =0. It is seen, as already pointed out by ‘Chwalla, that the introduction of 
primary bending moments results in a very small reductionof the critical load. 
it would, , of course, be hasty to predict that the relative insignificance of this — 
reduction applies to all problems of practical interest. In particular, when the ; 
‘geiamey bending ensaumeie are such that the el: elastic limit. of the material is ex- 
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ceeded, the effect of the | primary nary bending 1 moments ‘may y become far more pro- _ 
nounced (This is clearly true in the present example unless the members 
are assumed to be unrealistically slender. ) Conversely, it is not to be ruled ; 
3 that primary bending moments may actually increase the stability of the 
_ structure. This is demonstrated in the next section; attention is also called 
tothe well- -known fact that the lateral stability of beams" is increased under 
certain: circumstances if their ‘prebuckling deformation is taken into consid- 


MOMENT DISTRIBUTION 
~The form of the equations developed in the previous section suggests that 
a ‘moment distribution technique may be developed that permits the establish- 
ment of the solution through successive approximating steps. (Similarly, 
solution based on compatibility conditions corresponds to an 
distribution” method. No advantage is "apparent in the application of sucha 
method to the present problem; moreover, its convergence is less rapid and, 
in some cases, open to question altogether. ) For the simplified problem (that | 
is, in the absence of bending moments) this was first introduced by 
_ G, Winter, F. ASCE,13 P. T. Hsu, B. Koo, for the single-story frame and was 7 
extended by E. F. Masur, M. ASCE, 14 to multiple stories. In Winter’s paper!3- 
the magnitude of the load is assumed. d and the horizontal force necessary t to pro- a 
duce a sidesway is computed by means of moment distribution. So long as this 
force is positive, that is, in the direction of the assumed sidesway, the frame 
is stable; by repeated trials the critical load, which corresponds to sidesway 
without lateral force, can be obtained through interpolation. In Masur’s work il 
_ the process is generalized to investigate the positive definiteness of a matrix. — 
Again, interpolative techniques are employed, together with the establishment | 
of upper and lower bounds to aid in the solution of the problem. Terr 
_ In the present paper, the technique introduced by Winter!3 is extended tothe 
case in which primary bending moments are present. _ Further extensions are a 
~ possible, and fairly obvious, , for the multi-story building, to which the bound- 
edness principles of Masur!4 | can be shown to | apply also. A rectangular frame 
ike the one shown in Fig. 1 will be treated first; in that case it is permissible 
_to ignore the effect of the horizontal reaction H on the stiffness of the e beam 7 
; BC, although | such an effect could readily be included in the treatment. 
To start with, the magnitude of the force Q is assumed and the associated | 
‘rotations, -moments, and so forth, are determined by standard methods, as re- 
viewed previously. It is convenient to restate Eqs. 2, that, with the 


an 
K 


- 3 “Buckling of Trusses and Rigid Frames,” by G. Winter, P. T. Hsu, B. Kao, and 
-M. H. Loh, Cornell Univ. Experiment Sta., Bulletin No. 1948. 
“a 14 “On the Lateral Stability of Multi- ‘Story Bents,” by E. F. ‘Masur, Proceedings 
Separate No. 672, Vol. 81,1955. 
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in which the are b and cr refer to beams and columns, respectively, a 
and in which the effect of symmetry (that is p = 0) has been included. foo 

_ If now the structure undergoes a virtual sidesway Aw, the joint rotations 6, 


stiffness factor K, and the -over factor their value values slightly. 


+5 AE 
AM; +CA0;)+ 0; AK+ 0; (C K) - 
= which AK A K) are obtainable Eqs. 13, 14, 18. 


er 


in which the positive and negative signs refer to the right left column, 
_ When the additional ‘moments from the first of Eqs. 24 and from Eqs. 26 
are substituted in the equations of equilibrium (Eq. 11), the resulting system ; 
of equations is similar to the conventional slope-deflection equations, with the | 
term A Mf representing a fixed-end moment applied at the top and bottom of : 
7 each column, It is therefore possible to solve this system through a sequence | 
of successive approximations. In particular, the moment distribution techni- - 
_ que is most useful because it aneidie the determination of the additional mo- 
ments (and hence, of the additional shears) without the actual determination — y 
of the additional joint rotations. Owing to the linearity ‘of the system, the r 
sults depend on the magnitude of Aw. The therefore 


is 


a a value of Q and find all ‘stiffness factors and bending moments 
bat (2) _ Assumea sidesway of arbitrary magnitude, say, Aw = 1. This is — 
alent to the application of a fixed-end bending moment given by Eqs. 26 to the 
columns at each joint. It is noted that, in the absence of primary bending mo- q 
ments, this reduces to the procedure suggested by Winter. 13 It is also noted — 
-that no fixed-end bending moments are applied to the beams in nthe poesest case. 
. (3) - Balance the joints and determine the final _ bending moments A A Mf; int the 
columns, using, of course, stiffness factors that a are = by the presence - 
of the assumed axial forces, 
Compute the additional column shears rs from 


=A 


— 


— 
it follows from Eq. 24(b 
y 
— 
4 
q 
| 
: 


(5) Check the horizontal Loree AF 


AF =- 


b) If joint is fixed, then 


- ,2Q Q K’ Mo aw - Bin 
ij L in 
the present example, the initial of was assumed be equal 
or, 1. 769 E I/h2; then, with a = 1/3, » Mbc = - Mbp sind 
(0,222 @ L = - 0,393 E I/h for h = K, = 2. 630 E 
of 9, and in viewof aad symmetry of the initial deformation, the wisn bend- 
obtained by means of moment distribution otherwise). From Chang’ 
table? and by E Eqs. 29(a) and 29(b) it follows tl that at 
4Mp, 4Mop = [ (0.300) (0,223) - 2, 630] = raw; 


after. moment distribution and in consideration of the antisymmetry of the ne de- 


From Eq. 27, the additional sidieaian shears are 
"By Eq. 28, the force F is positive; 1 the structure is, therefore, stable and the 


assumed value of Q is too low. 


— error estimate may now be obtained peininatbeeiitie by setting F = 0 and 
iby solving a 27 and 28 for a using tt the v values of the additional bending mo- 


= tive, negative, or zero; this implies stable, unstable, or neutral equilib- a 
” respectively. The summation in Eq. 28 extends over all the columns in _ We 
special cases call for a slight modification of the procedure: 
= 
| 
| 
— 
— 
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“ments that were found previously. The value of Q= = 1.784 E I/h2 so so ‘deter 


mined, as proved for the simplified case presented by Masur, 14 represents an 
bound to the actual critical load; 


The next trial Q therefore lie between these two values. 


- 
As another example of the effect of primary bending moments the gabled 
frame shown in Fig. 3 is considered; for the sake of simplicity all members’ 


are assumed to have the same E I values. Again the moment distribution tech- 
nique described earlier in this section 1 is few modifications 


To the ettect of the position of the loads their value 


- the problem is first analyzed under the assumption that the loads are applied 
at the eaves as shown in dashed lines. s. In that case, no primary bending mo- 
ments appear and the type of analysis p proposed by Winter*” can be carried ov out 
without difficulty. Because of the symmetry of the structure e the problem - 7 
again divid divided | into possible modes of 
a that retains symmetry, 2 and 


(2) an antisymmetric involving sidesway of the e roof struc-_ 
= 


In the first case the critical value of @ can be shown to be determined by - 


\ KL 


&§ — 
— 
Ex 
7 
— 


in re o is given by Eqs. 3 with reference to either column and K is the as- 
Felationship: 


It is « easy ‘tosee see that thea antisymmetric c mode is associated with asmaller val- _ 
ue ¢ of and is therefore critical; the is of y and assumes the 

As the second approximation to the actual problem, the force is yone-ond ol 


= moments is to be. ignored. This may be visualized by i imagining a an ninextensible : 
— cable attached to the structure as shown in the figure. In that case the 7 
structure is under purely axial stress before buckling if the effect of the echange 


in length of the members is ‘ignored, 


_ Again two modes of instability are possible, ‘of which t the one associated with | 7 
_ the smaller buckling load involves sidesway. It is governed by the relation-_ 


% 


in which h& ¢ and Ko. refer to the column and the rafter, respectively, This le leads 


to the critical load Q=1. 10 | VL, that a substantial reducts tion 
the cross cable used previously is removes the problem is , drastically 
modified through the development of prebuckling deformations and bending mo- 
“ments. The analysis is lengthy and is therefore not reproduced here for a 
7 sake of conserving space. It follows roughly the outline of the analysis of the — 
rectangular frame, the chief modification being the inclusion of the initial ro- 7 
_ tations of the columns as well as the somewhat more complicated _ Kinematics 
(that is, two degrees of freedom). Also, the change in the axial force magni- 

tude in the rafters during buckling must be included. The results of the anal- 
ysis indicate that the structure has been stiffened considerably as a result of | 
the” prebuckling deformations. In fact, the first trial solution indicates that 

the critical load lies between 1.25 E 1/2 and 1.38 
sit is interesting to note me this seems to violate a principle of Rayleigh, — : 


apparent than real; two different neutral equilibrium configurations are being 
= here, and the problem falls therefore outside the pale of the Ray- 
leigh principle. (A similar “contradiction” exists in that a shallow spherical 
capor a shallow arch can be stiffened against snap-through by relaxing its 
‘boundary conditions i in such a Way | that it may undergo certain types of sym- 


_mations the load in the rahers ; is carried only partially in direct compression, 
‘the other part being taken u up in shear, . (There is at least some | question as to — 
7 "whether these results are entirely reliable. The assumed inextensibility of the 7 
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members an whose. should be 
_ investigated further. Certainly some obvious contradictions may arise. For | 
- example, if the structure shown in Fig. 1 were loaded by a concentrated force 
at the center of the cross beam the analysis presented here would apply with 
= small modifications. If now, in addition, a column were added to the structure 
a, directly under the load, then (in view of the assumed inextensibility of —_ 


this finite increase would take place for an arbitrarily weak center 
result that is patently absurd.) However, as in the previous example elastic 
buckling is highly unlikely for conventional slenderness aa Me 
CONCLUSIONS 
eo"  % appears that, with th the techniques introduced in this paper (and especially _ 
by with the moment distribution procedure described in the previous section), | ll 
- labor involved in determining the elastic buckling load of the actual problem 
has not been increased significantly over that involved in the “simplified” prob- 
dem. It is, therefore, open to question whether the customary neglect of the 


_ effect of primary bending moments is, in fact, justifiable from the point o of 
view of economy of computing labor, 


It is true that, for the examples considered, primary bending moments 3 af- 
fect the stability of the structure only very little; whether this is always true, _ 
however, is not clear. In any event, the stability of a partially plastic struc- 
ture is certainto be intimately related tothe presence of primary bending mo- 


ments. . It is ‘suggested that the tools described here may | ‘be useful in evalu. a 
ating this important problem. — — 2 
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Proceedings of the American Society of Civil Engineers 


ANALYSIS OF FRAMES LOADED INTO THE PLASTIC RANGE ae 


Aprecise method of analyzing plane frames stressed into the inelastic range 
_ is presented. The method involves the determination of a “compatible” ——— 
7 and rotation at a joint by the intersection of two moment versus end- rotation — 
curves s of the adjoining ‘members. The effects of initial residual stresses and 
: instability of beam-columns are considered in the analysis. It is shown how 
the behavior of frames may be predicted after the attainment of the ultimate 
_loads. Throughout the paper the frames are assumed to fail by excessive bend- | 
_ ing in the plane of loading. A complete example is included to illustrate the 


technique 


Two assumptions that are commonly made in the analysia of 


and that reduce the problem to a linear one in structural analysis, are that (a) 


the material mmees Hooke’ s Law and (b) the deformations of the loaded struc~ 


cited to Beene the formulation of the equilibrium equations on the basis of the 


_ Note.—Discussion open until January 1, 1962. To extend the closing date one month, - 
a written request must be filed with the Executive Secretary, ASCE. This paper is part — 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the — 
American Society of Civil Engineers, Vol. 87-8, No. EM 4, August, 1961. _ : 
Assoc. Prof., Dept. of Civ Engrg., Ohio State Univ., Columbus, Ohio. 
2 Research Assoc., Fritz Engrg. | Univ., Bethlehem, Pa. 
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publications 

-_ no provisions for the effects of the deformations or ‘the effects of : 7 

axial loads onthe moment-curvature relationship of the members. Other stud- 

bs ies that consider axial load and deformation effects are limited to cases in 
which the stresses do not exceed the proportional limit in any part of the struc- _ 

a * Two papers), 6 fall into this category of elastic analysis. They are pri- 


marily concerned with the determination of the sidesway buckling load of sym- 
metrical frames when the loads are not applied directly tothe joints. — 
At third method of analysis, 7,8,9 takes into account both inelastic action and 
| deformations, but in a somewhat restricted way. it ignores the effect of axial 
load on the moment-curvature relationship, and although the displacements of 


joints the of conc load are considered, the 


Ultimate. Load - Simple Plostc 


Buckling 
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Vol. 186, No. EM 3, June, 1960, p. 1. 
8*The Analysis of Partially Plastic Redundant Steel” Frames,” ” by 
Australian Journal of Applied Science, Vol. 7, No. 1, 1956, p. 10. Pe aa i 
*“The Elasto-Plastic Analysis of Two Experimental Portal Frames,” ” by J. 
Roderick, The Structural En: inser, ‘Vol. 38, No. 8, August, 1960, 245. 
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In spite of the good correlation with experimental results that other methods 
"provide, 3,4,5,8,9 the interest in the buckling of frames loaded into the inelastic 
_ range has made it desirable to develop a more precise method. The present 
method will consider inelastic action, axial load effects, and the exact ve seen 


the loaded structure. 


of specified cross sections is subjected to symmetrical loading as shown. Fer Pry 
a specified beam length and varying column heights, the curves give the ulti- 
mate load according to the different modes of failure assumed and according © 
te the methods used in the computation of the curves. _ (In Fig. 1 the load Pis © 
- equal to half of the total applied loads. It is equivalent to the axial force | 


_ a At the specified column height on the chart of Fig. 1, 1 line ne (a) gives the car- 
rying capacity according to simple plastic theory. ‘This theory assumes failure _ 
_ by symmetrical bending without sidesway. It ignores the reduction inthe mo- _ 
_ ment capacity of the columns due to axial load and due to the secondary os 
ments in the columns resulting f: from their deformations. 
Curve (b) assumes purely elastic behavior, with failure caused by sidesway 
~ buckling a as shown in the right sketch of Fig. 1. Points o on this curve are de- 


“ia by the methods presented elsewhere. 5,6 a 


paper. "Failure occurs by excessive symmetrical bending deformations; this 
mode of failure would occur when a lateral restraint at an upper joint prevents. 
sidesway motion. _ As will be indicated in an example, ‘the ultimate load = the 
largest for which an equilibrium configuration of the frame can be found. a 
_ Curve (d) gives the true ultimate load of the frame when it is not laterally 
' restrained. A procedure» is developed + elsewhere10 for determining a point on 
The method may be summarized as follows: 


mc. Assume a value for the load. te 
Perform. an exact analysis” according to the method described in the 
2 i) paper. The analysis will give the bending stiffnesses of the frame 


-members for the assumed loading, 
1 Perform a stability analysisto determine whether positive work is re- _ 
> "quired todeform the frame into an adjacent anti- -symmetrical shape. This may 
be performed by the methods outlined elsewhere! 1 (where only elastic frames 
considered), 
4) If step (3) indicates a higher or lower an load than the one as- 
lead in step (1), a lower or higher load, respectively, is assumed until the 1e 


_ “analysis of step (3) and the assumed load in step (1) are identical. Py wer 
Notation. _—The letter symbols. adopted for | use in this p sega are defined and 


_ arranged alphabetically, for convenience of reference, in the Appendix, © 


1 “Stability of Elastic and Partially Plastic Frames,” by L. W. Lu, thesis presented — 
to Lehigh Univ. in Bethlehem, Pa., in ‘ween in partial fulfillment of f the requirements for 
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The theoretical analysis is subject to the following limitations: pene enme 


The frame is symmetrical. 


4 2) The loading is symmetrical. 


(3) The members are prismatic. = 
(4) Only deformations due to flexure are considered. 
: o) The members of the structure are made up of columns that have hud 
loads applied to them except possibly at their ends, and beams that have neg- 
ligibly small axial loads in them. Although allowance for the axial loads in the 
_ beams may be made by a trial procedure, the present analysis ignores these 
effects because of their small influence | on the structural behavior. cong 
_ Of axial load in the columns is more critical and it is, therefore, taken into 
_ (6) The method is non-historic. It is necessary to specify that during the 
~ loading: there is no strain reversal of material stressed beyond the elastic 


a A typical structure, shown in Fig. 2, will be used as an illustration of the 
at method. The frame is symmetrical about a vertical axis teed 


= 


"through point C. It is desirable at first the columns and ths 


om In Fig. 3(a a) the column element BA is shown with axial load P and moment 
MBA applied toend B. Fig. 3(b) shows how the rotation of end B varies as_ 
is slowly increased while the axial load P is maintained constant. A pro- 
cedure > by which the curve in ‘Fig. 3(b) may be obtained» follows from t! the fact 
a that a column such as BA 1 may be considered as as a segment of a column deflec- : 
7 tion curve. The column deflection curve is defined as the shape that a com- 
‘pressed member (identical te to o the column member except in length) will take 
when held in a bent configuration by axial loads applied to the ends. . Such - 
_ column deflection curve is shown in Fig. 4. It may be identified by the angle a 
_ @q of the initial tangent. An infinite number of column deflection curves are 
possible fora single cross section and axial load. For each value of 8o, the 
column deflection curve will have a different shape. if 8 is small so that the 7 
stresses in the column deflection curve do not exceed the proportional limit, 
the curve is a sine wave and its half baa length is is determined by t the | Euler 


> 
= 
‘a 
le 
= 
fe 
(A: 
4 
4 
— 
4 b 
| 


‘section is equal to Py. It will ‘be noticed that the ‘segment AB of the column — 
fae curve in Fig. 5 resembles the column AB in Fig. 3(a). The only 
_ difference between the two curves is that the component t of the reaction at A in 

Fig. 3(a) along the direction AB is P, whereas in Fig. 5 it is P cos a ve 

» In — Visi the cia between the reactive forces at A and the direction AB 


It is seen that the angles between the total thrust at A and the direction AB are 
the same forthe segment of the column deflection curve ona the actual column. ae 

a Next the the magnitude of these thrusts will be examined. In Fig. 3(a) the thrust 


we. 
equal to 10% of 
_P, R exceeds P by 1/2%; and for smaller ratios of Vv to P, ‘R approaches P- 


very rapidly. ‘The preceding consideration makes it possible to consider the © - 


, a en of Fig. 5 equivalent to the column of Fig. 3(a) for all practical cases. 


on moment- curvature relationship of the ‘column is used here 

« - because it is easily adaptable to different ye sections and material proper- J 
ties. _ ‘The moment - curvature relationship i first for the > actual 
for determining the M-@ curve and performing the integration for the column 
deflection curve are outlined elsewhere. 12, For rolled steel wide-flange sec- 
tions with residual stresses due to uneven cooling, another “publication!3 pro- 
_vides information for the determination of the M-@ curve. _ The integration 
gives the slope and deflections at intervals along the column deflection core 
corresponding to the interval of length used for the integration. | It has been 
-noted!4 that with stresses below the proportional limit, the numerical integra- 
- tion gives lengths and maximum ordinates that are in error by less than 0.13%. 
_ In the comparisons given elsewhere!4 an interval of length equal to four times 


the : adius of gyration of the section was used. _ 
13% " Plastic Deformation ¢ of Wide- Flange B Beam Columns,” by R. L.. Ketter, E 
- Kaminsky, and L. S. Beedle, Transactions, ASCE, Vol. 120, 1955, p. 1028. Sle ae 
— Me Restrained Columns,” by M. Ojalvo, thesis presented to Lehigh Univ. in Bethie- 


chem, Pa., in 1960, in ees sat ae of the en for the degree of Doctor of i 
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PLASTIC RANGE 
After several column deflection « curves have been computed, the moment- 
rotation cu curve (M-@ curve) of Fig. 3(b) may be determined from them. Point 
Bis located on each column deflection curve at a horizontal distance h from 
one end of the curve (Fig. 5). The moment corresponding to MBA is Pyp 
while the angle 6 is the slope of the tangent at B plus yp/h. Each column 
- flection curve is thus seen to furnish the coordinates of a ‘point on the M-@— 
_ Beam member BD can now be considered for the purpose of determining an _ 
- Mpp -@ relation analogous to the Mpa -@ relation for the column. The con- ~ 
struction an MsBD ¢ curve for the beam ‘proceeds from the fact that 


4 
‘ 


center section of the beam is fully plastified (formation of a plastic hinge), - 
the Mpp -9 curve becomes simply a straight line parallel to the @ axis; Fig. 
10 for example. The deformed of the beam in 6 be referred 


Curve at t intersection 


to as a beam deflection curve. bending 
moment at Cc (Fig. 6) the bending moment at distance x measured from Cc is: 


A positive | bending moment is one that causes tensile stresses on the outside 

: of the frame. _ By using an integration ‘procedure similar to that used for the : 
columns it is possible to construct several | beam deflection curves, each fora . 
different value of the center moment ment Mc. beam deflection curve: 


value of 6 at x equal to is determined the integration. 
; Each beam | deflection curve, therefore, provides the coordinates of a point on — 
the beam curve - in Fig. 7. Because the axial load in the beam is 
small, a moment-curvature relation for zero axial load may be used for nu- 
_ merical integration without the introduction of an appreciable error. _A further ’ 
; refinement tothis part of the procedure would require that an estimate be made — : 
Vand that a slightly more appropriate M-@ curve be used for curve 
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When the two M- 9 curves, one for the column and one for the beam, are 

- plotted as in Fig. 7 the intersection will give the moment and rotation at B(or 
_ D) for the equilibrium configuration of the structure. In other words, of all _ 
_ the combinations of moment and rotation at B that are possible in the column — 

and of all the combinations that are possible for the beam, the correct one is 

the one that is simultaneously compatible for the beam and the column. Once 
_ the correct M and @ at Bare determined, it is possible to obtain the forces 
: in the rest of the frame by statics. _ Deflections may be determined by inter- 


-_ the basis of the Ragneses concept - of inelastic buckling, in which no unload- 
ing is assumed to take ee at the instant of buckling, the stiffness of the col-— 


umn is given by the slope of the tangent to the column M- 8 curve at its inter- 

section with the beam M-6 curve (Fig. 7). The stiffness Kpp or Kpp) for 
the member BD may be determined in the usual manner by considering the ef- 
“fective flexural rigidity (reduced due to yielding) of aa the ‘sections. s. Detailed 


= procedure for computation can be found elsewhere. 10 7 


thet behavior of the on frame shown in Fig. 8 will be investigated in ‘this i 


‘using the the procedure d described above. The dimensions and | member 
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1000;— Mp =950 ft.- 


> 
Fi. V vs. CURVES OF THE ‘COLUMN 
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we 3. 05 kips 


FIG. 10. _Mp BC VS. CURVES OF THE BEAM 
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“FIG. vs. JOINT RELATIONSHIP 


Ultimate w= 4. 
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aaa 12. .—LOAD VS. JOINT ROTATION RELATIONSHIP © 
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PLASTIC RANGE 
sizes of the frame. were chosen on arbitrarily. The : stress- strain diagram for - 
7 the e rolled steel members is assumed to he ideally elastic-plastic with a yield — 
stress level at 33,000 psi. The M-¢ curves that were used inthe computations | = 
were taken from another source.13 These curves assumed a well defined re- 4 
‘sidual stress pattern in which a maximum residual stress of about 10,000 psi 
3 occurred at the tips of the flanges. Fig. 9 shows Mpa -6 curves obtained else- 
-where.14 The Mpp -? curves for the beam ( Fig. 10) were taken from | — 
The : ratio of the force F applied directly tothe column top tothe distributed 
' load w 1/2 on the | beam, y, is arbitrarily maintained at 2.5. Solutions for sev-— 
eral values of w, corresponding to several il stages in the proportional loading © 
we The values of w that are considered are w * = 1,22, 2.03, 3.05 and 4.06 kips 
ft. The corresponding average compressive stresses in the columns are 
 P/A = 3960, 6600, 9900,and 13,200 psi. (The axial forces in the columns, cor-_ 
_ responding to these values of stresses, are 0.12, 0.2, 0.3 and 0.4 Py, in which - _ 
Py is the axial yield load (Py = A Cy) of the column ‘member.) By combining | 
: the curves of Figs. 9 and 10, asa of 6 and MB are obtained for each value 
will be noticed from Fig. that intersections of beam and ‘column M- 
‘curves occur in pairs until the load reaches a certain maximum value for the 
frame. At this maximum load the beam and column M- curves become tan- 
gent. Fora frame that is restrained against a premature sidesway buckling _ 
_ failure, a higher load would produce no intersection and thus indicate that the 7 
capacity of the frame has" been exhausted. ‘For | the frame of ‘Fig. 8, the maxi- 
mum value of w is estimated to be 4.11 kips per ft. . The maximem value « of wo 
- according to simple plastic theory is 4.94 kips per ft, indicating a reduction of 
8% of the load- ~carrying capacity due to beam-column action. ‘Fig. 12, which 
; _ shows the rotation of the joint at B with increasing load, gives a ‘graphic indi- 
cation of the non-linearity of the structural action,  ——— 
. 3 It is of interest to know whether the ability of the structure to support load» 
diminishes suddenly after the load attains its maximum value or whether the : 
structure’ s load- carrying» capacity is maintained until the deformations 
-comevery large. This information is obtainable from the second intersections 
of the M- 8 curves which trace the load deformation characteristics after the — 
J maximum load has been reached. Fig. 11 indicates that the unloading will ‘be 
gradual as long as the M-@ curve der the column has a long flat portion at or 
near the moment. ‘It is known | from: a study of M-@ curves for rolled 
remains below 60 and the average compressive stress is less than 40% of 
yield stress of the steel (it is assumed that the column does does 
except by excessive bending in its plane). 14 ; 
A new “method of elasto- -plastic analysis for rigid frames has been 
in this paper. It is based on the graphical determination of a moment and ro- 7 


_ tation at a joint by the intersection of two moment- -rotation curves of the ad- — 


joining members. In constructing these moment- rotation curves, the effects 
of initial 1 residual stresses and the instability of beam- columns —_ be included 
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™ more precise than those obtained by other currently available methods. © ; 
The chief importance o of the method is that it provides information by v which 
< problem of sidesway buckling of frames in the plastic range may be inves- : 
tigated. It can alsobe used to determine the carrying capacity of a frame sub- 
to axial forces in the columns, but restrained against sidesway 


Me 


pte 


= 


"movement. 7 It is also of interest when it is s desirable | to study, the nahenen « of | 


frame after the attainment of the ultimate load. 
The method described herein has been illustrated by reference toa 


— 
18 
me 
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= j. metrical frames shown in Fig. 13. In each case appropriate sets of moment- _ Reef 
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ments are then determined by the intersection of these curves. ah : _ 


os ‘Their Components” currently (1961) being carried out at the Fritz ae 
neering Laboratory of the Cvili Engineering Department of Lehigh University — 
under the general direction - of Lynn S. Beedle, F. ASCE. The investigation is _ 
sponsored jointly by the Welding Research Council and the Department of the 
Navy, with funds furnished by the American Institute of Steel Construction, 
_ American Iron and Steel Institute, Office of ‘Naval Research, Bureau of Ships, 
_ The writers express special thanks to ‘Theordore Vv. Galambos, A. M. ASCE, 
ae George C. Driscoll, Jr., M. ASCE, for their constructive suggestions that 


have been incorporated in this paper. a. 


E _-modulus of elasticity (kips per sq in. ); 
=~ = concentrated load applied to the top of a column (kips) “a 
moment of inertia about axis of bending | in. 4), 
ft-kips\ 
stiffness of frame member 
radian}? 


= halts wave e length of a column deflection c curve ve (in. 


pe — Mx = internal bending moments (ft- -kips); 


= joint moments (ft- kips); 


= axial load in (kips) ; 


resultant thrust in column (xips): 

V shear in| column, horizontal reactive f force at frame support 


= "distributed load on the beam (Kips /tt) 


= deflection at B ona ‘column deflection curve ( (in. 


ne 
7 
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= “maximum deflection of a column curve (in.) ); 
rotation of a joint (radians); 
= column deflection curve 


of a member (radians /in. 
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_ Proceedings of the American Society of Civil Engineers _ 


ANALYSIS OF NONCOPLANAR ELASTIC RINGS. 
‘By D. H. innate F, ASCE, W. C. Liessner, ? and Yen-Ming Yeh3 


_ The six scalar equations required to determine the components of the in- — 
_ determinate moments and forces in loaded or distorted elastic rings are de- - 


using algebra. Simple sign conventions and easily remembered 


analogy o: or to. ‘the shear and torsion analogy for the e analysis of co- 
"planar rings, or to a corresponding analogy using a single elastic centroid in i 
_ Space for noncoplanar rings. _ The analysis of a superelevated helical beam is 


_ illustrated, using cylindrical rather than cartesian coordinates, ‘The expres- ee 
_ sions are easily converted to any orthogonal coordinate system. | 


ay ot “ 


The analysis of irregularly shaped elastic rings stressed by either applied 
we or by distortions is based on either the classical analysis of least a . 
or on the more modern concepts based on satisfying statics and geometry. By ; 
using the latter method Hardy Cross, F. ASCE, developed the column analogy4 4 
on - Note.—Discussion open until January 1, 1962. To extend the closing date one month, © 
8 written request must be filed with the Executive Secretary, ASCE. This paper is ml 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the — 
3 American Society of Civil Engineers, Vol. 87, No. EM 4, August, 1961. = = ne 
_ 1 Prof. and Head, Engrg. Mechanics Dept., Virginia | Polytechnic Inst., Blacksburg, Va. 
2 Engr., Sperry Gyroscope Co., Great Neck, N. Y.; formerly Instr., — Mechan- . 
ics Dept., Virginia Polytechnic Inst. Blacksburg, Va. | 
a Engr., Amman and Whitney, Cons . Engrs. » New York, N. Y.; formerly Graduate _ 
Asst Engrg. Mechanics Dept., Virginia Polytechnic Inst., Blacksburg, 
_ 4 *The Column Analogy,” by Hardy Cross, Bulletin 215, Univ. hed Mllinois Engrg. Lal 
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for coplanar rings stressed by forces or in the same 
plane, or by ‘moments or rotations acting about axes perpendicular thereto. 
fae P. Michalos, F. ASCE and Frank Baron, F. ASCE extended the method - 
in their shear and torsion analogy® to coplanar rings stressed by moments or 
rotations acting about axes lying inthe same plane, or by forces and displace- 
= ments acting perpendicular thereto. They also extended the method to cover 
elastic rings curved in space. 5 The special case of a helical beam without 
_ superelevation was analyzed by Alan C, Holmes, 6 Inall threee of these studies 
the analysis was based on a scalar rather than ona vector or approach. Other» 
_ studies on more general noncoplanar rings, such as pipe lines curved or seg- 
mental in space, have been made by John E. Brock’,8 using matrix analysis. 
7 aa Equations based on a vector or matrix derivation must, of course, be ex- 
- _ panded to the more familiar scalar versions when the actual computations are . 7 
_ made, Nevertheless, the vector derivation presented herein is considerably © 
i shorter and ‘the standard vector § sign ‘convention easier — to follow than that 
usually adopted in scalar derivations. Furthermore, the equations in vector 
_ form are more easily expressed in the coordinate system most suitable for a - 
particular ring. For example, cylindrical coordinates would usually be 
- appropriate for a helical beam. Other than this, however, the final resulting — 
scalar equations do not differ from the more familiar scalar analyses4;, 
_ those derived herein do present expressions for a single elastic centroid or 
_ The elastic axis of the ring considered herein, can have any arbitrary _ 
_ shape in space, and the modulus of elasticity and cross section can vary in - 
manner along the ring. The analysis herein is restricted to the case 
: which (a) the effect of deflections due to longitudinal strain can be neglected, 


; or superimposed, (b) the deflections are small, (c) the curvature per depth 7 


trary m manner > by either forces and moments, distorted at any ny point by relative 
displacements or small rotations, or throughout by temperature 


: _ Consider any continuous, ‘noncoplanar r g, subject tothe restrictions pre- 

" viously listed, in equilibrium and loaded i in any arbitrary manner, If the struc-| 

: ture is made statically determinate by cutting it at some point A (Fig. 1), the 
_ ends and tangents to the elastic axis at A will undergo’ a relative linear and 
_ angular displacement due to the applied forces andmoments. These loads will 


produce a statical force and moment which act on the cross section of 


% and moments M; at the ends Cot rigid brackets connected to the cut ends of 


5 “Laterally Loaded Plane Structures and Structures Curved in Space,” by James P. 
Michalos and Frank Baron, Transactions, ASCE, Vol. 117, 1952, p. 279.5 
3 ’ “Analysis of Helical Beams Under Symmetrical Loading, ” by Alan C. Holmes, Pro- 
ceedings, ASCE, Vol. 83, No.ST6, November, 1957. 
_ 7 “Matrix Analysis of Flexible Elements,” by John E Brock ceedings, First U.S. 
“Natl. Congress of Applied Mechanics, 1952, p. 285. 
8 “A Matrix Method for Flexibility Analysis of Piping Systems,” ” by John E. Brock, 


Transactions, ASME, Journal of Applied Mechanics, Vol. 74, 1952, p. 501. 


— 
— 
— 
( 
— 
tans 
— 


the 1 ring at A. . The heostete are coincident at no load. The origin n of the coordi- 
nate system is located in any arbitrary manner at end C of these brackets. _ 
+The restoring reactions will produce a restoring force Vr and ‘moments My on a 
the cross section of each elemental length ds of the ring. Hence the final mo- . 
ments and forces acting on a cross section of ds | 
M = Mg - i 


and the of both are satisfied, 
If, instead of being subjected to actual loads, the ring is distorted by te aes 

jecting the two cut ends at point A) in Fig. 1 to a linear displacement, or - 
subjecting the tangents to the elastic ; axis, there to a small relative angular 


1.— LOADED RING WITH REDUNDANTS- 
APPLIED AT THE ORIGIN TO RE- 
STORE CONTINUITY 


applying a suitable force Ve and moment Me at 


_ Sign Convention.—The rectangular components ‘parallel to the asain 


ends” and tangents at A can be displaced in any prescribed 
3 axes of the moments °¥, forces V, rotations dé, and displacements 6 will be 


a. considered positive if the vectors which represent them act in the direction of 


- the positive axes, provided such components act on the front face of a cross 


section of the ring. Normally rotations do not obey the law of operation of a 
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q 

a 

— 


vector even though they possess magnitude and direction. rotations 
: 7 g infinitesimal they can bet treated as vectors. The right hand rule is used to 


‘The front face of an element is defined as the cross section normal to the 
elastic axis farthest counter around the ring when viewed from the 


Face + 


thie: Retiating Linear ond 
Mone 


~ 


- positive end of a coordinate axis. Fig. 2 illustrates the positive directions of 

_ the vector quantities and positions of the front and rear faces. he ss 7 
_ Direction Cosines. —The orientation of any differential length | ds along the : 

elastic axi axis: x] or of the principal ax axes 8 X2 al and x3 ¢ of the ‘rings cross section 


as 
— 
| 
ii 
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. - perpendicular to the elastic a axis, canbe specified most conveniently by me means 
_ of direction cosines. Fig. 3 indicates the relative orientation of these axes as 
well as the unit vectors in the two rectangular coordinate systems, The axes 
associated with the element will herein be considered as the primed x’, 
and the arbitrarily chosen axes x, y, and z as the unprimed set. 
_ For any vector, such as one for moment, its three orthogonal components 
‘in the primed coordinate system may be related to its three orthogonal c com- — 
; ponents in in 1 the e unprimed system in index rotation Mt pe ath 
Mj = by Mj i, j= x, y, zor 1, 
in which Mj is the ‘moment vector M with. components in the primed coordi- | | 
nate system; M, denotes vector M with components in the 


angle between the x} and 3 x} axes. The vectors and direction cosines are ex- 
a in seed index notation. An unrepeated index (such as i given previously) 


as or, briefly, 2, repeated index (euch as as j) implies 
- summation over the range of that index (such as x, y, z). Hence, Eq. 3 can : 


can be written in transposed form as" 
’ 


“KY 


M,ds Mjds 


| 
| 
q 
- 7 __. The slender elemental ring segment ds, whcn subjected to bending moments a 
Ingle changes as follows 
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Herein, G E I'2, and so forth, to the s stiffness. The 
term J'is not a polar moment of inertia, except for for a ring of circular cross 
_ Elemental angle changes d6; may be related as was the moment vector Mi 


A jk M; ds ds = = by jk k1 Myds 


_ when Eq. 8 is expanded in x, y, and z, three scalar equations iaeaai to . 
following result, but only one is shown herein 


iM 


the: terms in brackets S are denoted by 2 and by etc. 


+ + 
xx 
Eq. 8 m may, then, written as 
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The e 1/Ky terms may be regarded as widths of elastic areas ds » tone. —_— re- 


lationship of these nine terms in terms of bij and Kj are listed in Eq. 41, 
Appendix I, It can be shown that Kij = Kjj_ and so only six elastic areas are - 
zs The expressions derived herein for dAjj differ considerably from — 
those proposed originally by Michalos and Baron,9 essentially their 
simpler functional relationships to the direction cosines. — 
 Statical and Correction Reactions.— The reactions Ve. = Me 
| (Medi, applied at the ends C of the rigid bracket of Fig. 1 which coincide with © 
- the origin of the coordinate system, will produce on a cross section of the ele- 
- ment ds, a restoring force Vr = i, ¢ and a moment which may be expressed in 


~< a4 


in index notation as 
in | which Cijk is a permutation symbol fora vector | cross product; Cijk = = +1 if i 
i, k can be permuted to 1, 2, 3 by an even number of interchanges; Cijk =] if : : 
the permutations are odd; and Cijk = if any 2indicesare the same. 
Eq. 12 or 13 may be expanded ir in scalar form and the enpeianse dropped 
for sake of simplicity to yield Eq. 43 of Appendix 
Corrective Angle Changes. ~The statical moment M, = and correction 
both produce differential changes (ae s)i' and (d at the 


is rei restored. Res. 11 and 1 


Tag \ 


ee now, the parenthesis for individual vectors are dropped for sake of sim- 
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rigid brackets at C due to the statical moment (Mg); and will be considered — 
as a load Pi on the elastic areas Ajj. The equations for P;, listed as the first 
‘three Eqs. 44 of Appendix * and id the the Ajj, Cij terms are defined in Tables 1(a), 


_ The left hand side of Eq. 15 ae a relative angle | change between the 


_ Corrective Displacements. —The differential angle Po a at the ele- 


and are substituted in Eq. 


é Expansion of Eqs. 19 yields the last three of Eqs. 44 of Appendix I. The in- 
dividual integrals ali Sify _ and Dy are defined in Tables 1(b), 1(c), and 


__-The left hand side of Eq. 19 represents a relative displacement of the ends 
of the rigid bracket at C due to the statical moments (Mg)m and will be con- _ 
sidered as a moment Dj onthe elastic areas 
i The Pj and Dj terms represent displacements at the origin C and may — 
- from actual distortions such as those accompanying shrinkage, support settle- | 

ment or rotation, temperature changes, etc., as well as applied 


| 


p.\. 
tions occur, such as those for 


in which a is the coefficient of thermal said 1; denotes the projected = 
_ length of the ring along the ith axis; and AT represents the Sempersmne change. 


FINAL EQUATIONS Al AND ANALOGIES 


The six equations for Pj and (Eqs. 4 44, Appendix represent relat on- 
“ships between these angle changes and displacements at C, , due to the statical 
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ELASTIC RINGS 
moment ( (Ms)i, in terms of the restoring forces (Vc); and moments (M,)i at Cc 
and the area elastic areas a and their first (Cjj) and second(Jjj)moments. 
Tf ce certain of the Cy terms are made to vanish by locating the origin of the 
unprimed system properly, the six equations may be simplified somewhat. | 
_ The equations then reduce to several popular analogies4,5 and become simpler 7 
forcoplanar structures; 
Column umn Aralogy.- —Consider a coplanar ring, symmetrical about and lying» 
in the xy plane and loaded by forces (or displacements) in that plane, or by 2 
moments (or r rotations) about the z axis. Let t the z axis coincide with the Zz 
axis, Then z = 0,bxz = = byz= bzy= 0, bez = 1, and byy = -byx. Since 
bending takes place about only the z' orz Ixis M, = = 0. Hence, allterms 
Eq. 44, Appendix I, vanish except those ‘solid lines. If the 


remaining three equations are further simplified by setting Cex = = Czy = 0, cor- 


Fecting the Jjj terms accordingly, and substituting solutions for Vey, 
a & Vxe in the last of f Eqs. 43 of Appendix I, the familiar column analogy 
q 


‘Shear and Torsion —Consider a coplanar | ring lying in and 
metrical about the xy plane and loaded with forces or displacements along 7 
the z axis, or by moments 0 or rotations the x-and y axes. Then, = 0, 
bxz = = = bzy = 0, bzz = _ and by as before, , but now M, =0. 
"Hence, all terms of the six equations vanish except those overlined with dotted 
“lines, Letting Cy, = Cyz = 0, the and terms accordingly, 


__ Noncoplanar Rings. —When Baron and Michalos developed their analysis for 
noncoplanar rings, they set all Cij except the Cij terms « equal to zero, In this” 
way they located three elastic centroids, one in each coordinate plane. These 
did not coincide with the projection of any single point in apace. They did 
_ Simplify the six equations so that each nowhad but four terms, but the P dij, ry 

- and D; had to be corrected for the transfer of the axes. 
similar simplification was developed? by setting the three | Ci =0 
corvasitnn the Jij, Pi, and Dj accordingly. This amounts physically to locating | ; 
a single elastic centroid or neutral point in space, and reducing the six equa- — : 
tions so that - each now has but five terms. Physically, this restriction means 
that a restoring moment Mex at C produces no displacement of C along the x 
axis, that a force V ex produces no rotation about the x axis, etc. | ae 
ae general, the labor of locating the origin at the centroid of the elastic 


areas, and of correcting the Jj; “moment of inertia” terms by analogous paral- 
lel axis etc., is seldom worth the effort, This 


ring is divided, assuming that the entire operation has not been soomnechy 


en on a computer. Tables 1 (a) to 1(e) illustrate the manner of deter-_ 


— 


9 “A General Procedure for Analysis of Elastic Rings in Space,” by Walter C. Liess- 
ner, thesis presented to Virginia Polytechnic Inst., at Blacksburg, in October, 1957, as —_ 


a 
2 
4 
2 
| found tO he identical with the shear and torsion analooy of Baron and Michalos 
7 
— 
2 
| 
i 
a } of any ring, coplanar or noncoplanar, is most readily _ Bs 
if 
— 
partial luiliiiment of the requirements lor ule degree Ol Master ol science. 
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ignored, because some subscripts are interchanged to conform to the index a 
notation used herein, and because the signs of some terms vary. These sign 
changes result from adherence herein to the right hand rule used in vector 
Tabular solutions in other than cartesian | coordinates follow “similar 
"patterns, The problem subsequently illustrated lends itself to cylindrical 
coordinates but a double transformation between three sets ortho- 


Note: Angle of -superelevation 


Angle o the 


« 
- NIT VECTOR ORIENTATIONS: FOR A ASUPER- 
ELEVATED HE LICAL BEAM 
ist broken into ai into at least ten segments. Tabular solutions are decidedly simpler 
than those involving integrals when the cross sections of the ring vary, or when 


i 


the elastic axis cannot be expressed in convenient functional form, 
_ Cylindrical Coordinates. —Consider a helical beam fixed at both ends and 
Let ajj = the direction cosines between (primed) cylindrical 


¥ 


i 
mining the coefficients of the si fersfrom 
that nrecented by Michalos and Barony hecause the centroidal corrections are 
— 
the summations involved in tabular are not as precise 
as integrals, the results are invariably quiteaccurate, especially whena ring 
‘a. 
| 
— 
= 
— 
tm 
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and bij direction cosines the (double primed) principal 
tional and bending elastic axes (1” ") a ) and the (primed) cylindrical 


- cos @ sin cos acos 
a - sina sing sin a cos 


= ay Dim =a, dA! 
The terms in Eq identical with those of 41, Appendix I, 
provide the expansion is now conducted with j in terms of a 1", 2", and dil "in-- 


stead of and and k int terms , and instead of x , and z 


,M 


pla em mas be 


mn 


3 
; 
— 
vector may also be expressed in terms of the components of the moment 
vector in the unprimed (x, y, z) system as 
dd. Af, b,. a... M_— ds = a,, ine 
(27) 
= f@ (28) — 
4 
= 
Hence 
Similarly # 
— = 
6 mF = fe. Dp: (ae fe. p. a,, dA 


= 


_ worthwhile to rearrange these tables so that the Ajj, Cy, and Jjj coefficients © 
“ean be more came computed. Tables 2 and 3 define such terms and show 4 


= cosh 


=z cosh 


Ez ( 


= | 


how. they are combined to form the coefficients of the final matrix (Eqs. 44 


helical ‘beam is fixed at two ro ends ow apart and loaded 
an upward ‘unit vertical load at its midpoint, Determine the variation 
bending and twisting moments along its elastic axis axis if as = 5, 6 = 30°, h= = 10 ft, 


), 32, and 33 could be pany of the t 
4 metitionc nature 
ae 
99 | C309 | C409 
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c. 


y ‘104 * * 


110 208 209 410 — S04! 
Gos - 2C210 + r (C505 - C303) 


~ 
Px 


ry = C110 209 Caio (C305 + Cso3)_ 


— Cy 


C407 


»,* 41, eon I, the Kij terms can be computed. Letting the 
and z' = = x3 and using 1, 2, and 3 instead of x, 9, and x 


— ft, K] = K3 = 10° 3 = JU” K-It". (2) Determine the 
— 
— 
By 
= + + (C399 + Ccoig) “604 “211 ‘ 
stituted in 


Length F 
“tint 
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indices for both the double | primed and primed axe axes, ‘it is found, for ommete Poe 


a The beam into t ten equal segments keep each term constant 
for each segment. Tables 4 (a), show the required 


Values of (Degrees 
| 


90% 10. 4683 
/ IFront Face 


to Getermion the coefficients o of Table 2 and the terms of Table 3. These final ; | 
= Cij, and Jij coefficients, and the Py and Dj terms found in Table 4 (c) are 
shown in Table - 4 (e) as they would be substituted in Eqs. 44, —_— 


3 
} 
on 
— 
+ 
i 
J 
Mi(ftib) | 0.40 
0.230 
— 


1961 


substitution of these values in Tables 4(f) and (g) yields the restoring moments 

— (My);. | The final moments and shears are found by using Eqs. 1. Results 

moment are plotted in Fig. 5 in terms of M,/Fr and Mj /Fr versus ¢%, and the 
values of Vi and Mj are shown on the ends of the helical beam. (b) The stiffness 
factor about the y axis at A is defined a as the moment at A required to rotate | 
the helical beam about the y axis at A through a unit angle with all other sll 7 


Ay= 


Vex) cx 


4 


=1, , and so forth, may be obtained from —- after 
Eqs. 26 and 27. 
i 


9a) 


11. 03. 


+0. 03642883 M... 0801448 


= -0.03022609 M = +0. 3491825 
V_=-0, ).01457716 M = - 0.4658233 
May 0=+0. 1884 = "Say 


All values for the preceding Wodir and terms must be 
“4 by 109° since the 10-> factor was deleted from the 1 Kij terms as used in 


BSolution of th 2 yaticr vi a the rec te | 
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LUSIONS | 
The herein for the six equations needed to deter- 
mine _the redundants in a loaded or distorted statically indeterminate ring 
| in space, considerably facilitates their derivation and uses standard ~ 
mathematical sign conventions for moments and rotations. Modification of 
these equations, when expressed in index notation, is especially simple if, 
~ auxiliary coordinate systems need to be used in the analysis, ‘aati the case 
forhelicalbeams, 
_ It is shown that the equations reduce to the column analogy ont to the shear 
_ and torsion analogy for coplanar rings. However, for noncoplanar rings there 
_ seems to be little advantage in modifying the equations so that the redundants 
can be located at the elastic centroids. mrecisag 


“Analysis of Helical Bar Type Structures in Coordinates, by 
Yen-Ming Yeh, thesis presented to | Virginia Polytechnic Inst., at Blacks- 
burg, in May, 1958, as partial fulfilment for the degree 
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F EQUATIONS IN CARTESIAN COORD 
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XX XZ 


ay = 1 radian, the units are kips for the (Vc); terms and kip-ft 
"ADDITIONAL REFERENCE | 
b2 be | 
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+x 
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= cz XX CX “XY Vez 
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for distorted struct ure: 


m gusts 
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MOVING LOAD ONA A ‘LAYERED. HALF PLANE 


a surface nuclear burst is that of the uniform motion of a steady pressure dis- 2 
tribution over the surface of a half space. In this paper a special instance of 
that problem is considered. The steady state displacements and stresses pro-— 
ducedin a two-layered elastic half plane by the uniform superseismic motion 
of a constant, concentrated normal line load along the the surface of the half plane 
are determined interms s of dilatational and distortional potentials. The solution 
is constructed by systematically adding to the solution of the analogous prob- 
lem for the one- layered half lf space the solutions for multiple reflections and 
‘refractions. of plane elastic waves from plane boundaries. ‘Superposition of the 
_line-load solution can then be utilized to obtain the solution for an arbitrary 
steady pressure distribution uniformly progressing with superseismic speed 7 


Notation. —The letter symbols used in paper are defined 


“appear, in the illustrations or in the text, and are arranged ieninenaiin 


convenience of reference, in the Appendix. 
Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
-American Society of Civil Engineers, Vol. 87, No. EM 4, August, 1961. — 
‘Asst: Prof. Civ. _Engrg., Univ. of of California, ‘Calif. 


i 
Sackman,* A.M. ASCE 4 
— 
A nroblem of interest in studies of transmission of cround shock cauced by 
4 
if 
(oe for a particular medium subjected to a step pressure of uniform intensity and 2 
step pressures which decay exponentially, 


_ Angeat, 


half space the uniform of load of con- 


stant intensity alongthe surface of the half space are investigated. The speed © 
.- the load is assumed to be greater than the velocities of propagation of dila- 7 

tational and distortional waves in either layer (superseismic case). It is fur- | 
ther assumed that the load has been moving steadily for a long time sothata — 

steady state (of plane strain) exists with respect to a coordinate system at- — 


_ tached to the moving load. The solution resulting from a a steadily moving arbi- 


tion, ,and as a step with muameain’ decay. Although not considered herein, the 7 
case of tangential loads on the surface of the two-layered half space can be _ 


SOLUTION» 


AN 
LAYER 


has adepth h, and the bottom layer : an infinite: depth. The c concentrated vertical © 
line load of intensity P (measured positive downward) moves over the a 
of the half-space inthe negative x x direction with a speed U. At time t=! 0 an 


two coordinate systems are taken | to be in coincidence, so — a _ oa 

The solution will be given interms of dilatational (or longitudinal) and distor- a = 
1 tional (or transversal) potentials, 6 and W, from which the are 


— 
toa stripution Can be odtained DY Superposition OF the line 10aq Solution, 
Let (x, z) denote fixed spacé coordinates and let (x, z) denote space coor- 
dinates attached ta the moving Inad P (Rio 1) Theton laver of the medium 
— 
— 
— 
lig; | 


Strain and stress components are then derived from the displacements in the 
ual manner, Superscript (1) will identify quantities in the top layer, and su- 


rs script (2) quantities in the bottom layer. 


The solution is obtained by utilizing the results of the analogous problem ci 
for the one- layered half space? in conjunction with the results for the reflec- 
and refraction of plane elastic waves at plane boundaries.3 ‘Fig. 2 shows 


the situation for the one-layered half space, in which the solution consists of 
dilatational wave of strength anda distortional The 


us 


has the ‘simple character shown in ‘Fig. 2 for the 
case, that is if the Uis than velocities of wave 


-Toforma a two- -layered half sp space, a plane interface at depth Zz: =h is intro- 
— and the two layers are assumed to have different properties, Fig. 3. 
plan 


ge 


2 *g Plar o] 
n al of Applied Mechanics, Vol. 25 , No. 4, December, 1958, pp. 433-436. samen ina 
“The Earth, H H. Jeffreys, “Third Ed., Univ. Press, New York, 1952, 


— HALF PLANE 
, Inthe upper layer the 
| 
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ow 
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generated, A similar phenomenon occurs due to the distortional 
The waves signified by gly and vy “must, ‘in turn, be reflected off the | plane 

boundary at z = 0 so as not to effect tractions on this | boundary. 
flections occur -ad infinitum. Fig. 3 depicts the beginning of this yracene, The 

0 


notation of N. Davids4 is adapted for use in this problem, so that t ol, vi 

‘represents the strength wave inthe top lay. er gen- 
erated from iprevious dilatational waves andj previous waves after 


successive reflections. A A similar notation on » it+j = odd integer 


is for the of waves in the layer. 
laws of reflection and refractionare nowused to determine the strengths 
= of the generated waves in terms of the strengths of the waves incident onthe _ 
= 0, h. When incident waves impinge on the interface at 
= odd integer; “i 


Aca) | 
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ij“ 


alt) 


A=: 


+b 
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ip... 


— 
| 
4 
— 
Ad A= (a, +5) - (4 %) 
ob Applied Meshanis, Vol. 28, No. to te 
» VOL. 6, No. 4, D on in Plates,” by N 
» December, 1959, pp. 651-660. Davids, Journal 


Fig 


4 
, . , are the strengt 
note the strengths of the resulting reflected waves, re 
gre aves, U is the speed of the moving load, 
, C,, refer to the respective velocities of propagation of dilatationaland 
_ distortional waves inthe layers, odie , HW ~ are the Lamé constants of the lay- 


rs, andp is the mass density of the layers. = + 


reflections from the free surface = 0) the (recurrence) r relations 


are (valid for i + j= even — 


- (at. 
m 


pl?) a?) . 


at a matter of notational convenience, ote: recurrence relations Eq. ,3(a) 
and 4(a) may be unified as 


1 


-(t. an index is sagen, the cor rresponding a is zero. ) In this way an 

of the strengths | may finally be related t to the two primary strengths, 9 

4 *o » given by the | solution of the analogous one-layered problem: 


g 
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Adding the effects of all the reflections the solution to the 
 two- problem is finally 


= 


The is thefunetional form form of the primary ry generating way waves of 
yh) 


) given by the ‘solution of the analogous | one- layered 


F' is the derivative of F, and H is unit step “(it is ‘not. 
necessary to determine F, as only ’ derivatives of F are required for the dis- 7 
placements.) Interms of the potentials éand W of Eq. Ta), the displacements 


=2y 


4 
= 
ae 
— 
— 
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form | step 5 pressure, Po H(x), Fig. 4 (a) is given in | Fig. 4(b) for a a point | in a ; 

top layer (x = 0, z = h/2) and in Fig. 4(c) for a point in the bottom layer 

(x = 0, z = 2h). The material properties have been chosen such that % =u, 


(1) 

in which sound velocities ar are higher in the bottom layer than in the top layer, 

a situation likely to be encountered in nature. _ For this same ‘medium Figs. 5, 
6, and 7 depict the early time rene histories progaced by an exponentially 
decaying step pressure loading, py exp (- «x/h) H(x), in which pp isthe maxi- 
-mum value of the appliedsurface pressure andx is adecay factor. The results 


3 for three values of xk are illustrated: = 0. 1. 0, 5.0. 


CONCLUSIONS | 


‘ai stress distribution and displacements produced by the uniform, super- 
seismic n motion of a concentrated line load of constant intensity over the sur- 
face of atwo-layeredelastic half have been determined in terms of dila- 
tational and distortional potentials. The numerical examples considered treat 
a half space in which sound velocities are greater in the bottom layer than in a 
the top layer. For such a layered medium, and for a uniform step pressure _ 
progressing over the surface of the half space, the peak value of the stress, in 
either the top or bottom layer, occurs after several reflections and refractions — 
of the elastic waves generated by the moving surface load. These peak values © 
are greater thanthose which would occur in aone-layered half space composed 
of either the material of the top layer or of the bottom layer of the two- layered 
For the case of decaying step pressure distributions in uniform motion along 
the surface of the layered half plane, the peak stress at a point in the medium 
- occurs at the arrival of the first signal caused by the surface load. Thus, | for 
loads with high enough decay factors, the peak stress in the top layer can be 
_ determined fromthe (simpler) solution of the one-layered problem, 2 whilethe | 
peak stress in the bottom layer can be obtained from the consideration of only - 
a single elastic wave refraction at the interface of the two layers. In this case 
the maximum stress in the top layer is the same as that which would occur if 
the medium were one-layered and composed of the material of the top layer. 
The maximum stress in the bottom layer is again greater than those which 
would occur in a one-layered half space composedof eitherthe material of the — 
top layer or of the bottom layer. 
It is noted that as the decay factor of the loading increases, the interaction | 


at a point in the medium between the various reflected and refracted onatic 
waves” decreases, Hence, for most ‘Points in the medium, if the surface load 
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has a highdecay cay factor, by the time a a signal from a particular reflectedor re- 
fracted wave is received, previous signals wil will have decayed to a low level, 
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APPENDIX, —NOTATION 


a: The e following letter s 


= wave constant ‘defined Eq. 3(b); 


= constant defined by ‘Eq. 


by = constants ned by E Ea. 


velocity of dilatational ar and dit distortional wave, respectively, 
‘defined by . 3(d); 


= distance defined by 10); 


- a functional form of the primary generating waves defined 
by Eq. 8; a 


= 
= “unit step | function; 
= depth of top rs 


ration of shear moduli defined bs by Et Eq. 3(d); 3 


Mach number of "moving load with respect to dilatational 
and distortional wave velocities, respectively, defined ial 


~ of the Mach | defined by Eq. 


velocity of moving surface 


= = moving rectangular cartesian coordinates defined iby Eq. 


= fixed rectangular cartesian coordinates illustrated in Fig. 


= - angle of inclination of dilatational and distortional plane 
Fespectively, illustrated in Figs. 2 and 3; 


= constants defined by Eqs. and 4(b), 


5 

ia August, 1961 iM¢ 69 

— 
a r use in this paper: _ 
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H 
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= normal distance behind associated distortional and. dila- 
- tational plane wave front, respectively, defined by Eq. Tb); 
= Lamé constants of the saburtels 
= mass density of the material; _ 
= cartesian components of the stress tensor; ’ 


and nd distortional potential, respectively; 


Teflections ‘and/or refractions; 


Laplace’ s operator; 


= 


Qh 
in 
"ie 
q 
be 
“iy 

| 
superscript (2) “= quantities in the omiayer, 
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CONJUGATE BEAM SIMULATED BY ELECTRIC CIRCUITS | 


ByS. Freiberg) 


is made with electric circuits using strings of resistors i in which voltages, 
currents, and resistances simulate moments, forces, and lengths. The avail- — 
_ able voltages, simulating moments, of the “electrical beam” are used | to pro- 
—_ duce currents, now simulating forces proportional to moments, used again to 
_ load another esectrical beam of the same “length” as the first om, which will 
_ behave, therefore, as a a conjugate beam. _ Examples of circuits, consisting only 
of resistors, are given for typical cases together | with experimental results; 
_ deflections of beams, given by the conjugate beam, are determined by this 
method with — errors W when compared with deflections computed analyti- 


cally 


us 
- a In an earlier paper, 2 analogies between electrical resistance: networks and 
beams were analyzed. Forces, moments, and lengths on on the beam are simu-_ - 
lated by currents, voltages, and resistances, respectively. The elastic prop- — 
erties of the beam, necessary for the resolution of indeterminacies, are intro- 
“duced by an auxiliary which applies moment area theorems. 


Note.—Discussion open until January 1, 1962. To extend the closing date one a 
a , wedtion request must be filed with the Executive Secretary, ASCE. This paper is a 7 


of the copyrighted Journal of the | Engineering Mechanics Division, ee of the © 
American Society of Civil Engineers, Vol. 87, No. EM 4, , August, 1961. des = 


1 Prof. of Electrotechnology, Universidad Nacional de Tucuman, Argentina; ne 
Prof. at the Univ. of California Div. of Mechanics and Design, Berkeley, — 

rs 2 “Electric Circuits for the Simulation of Bending in Beams,” by S. Freiberg, Paper — 
presented at the S.E.S. A. spring meeting, 1960.0 
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The analogy ‘indicated above is also used and adds a complementary 
circuit for the simulation of the conjugate beam. 
This “electrical conjugate beam” can be realized in the same way as the 
auxiliary ‘beam used to apply the moment areatheorem, thus providing a use- 
ful tool for solving rapidly determinate and indeterminate beams, sisal 


moments shears, and by of voltages” and currents. 


The . well known conjugate beam is 5 defined as an a beam of the : 
‘ome length as the real beam and when loaded with the M/E I (moment divided — 
: by rigidity) diagram of the real beam, its new shears and moments will ‘repre- 
sent the slope and deflection of the actual beam. = 


It is known from statics that in a any b beam —_ 


in which y is the deflection of the elastic curve; 6 denotes the slope of the | 
elastic curve; ; x is the abscissa on the beam; M represents the bending mo- 
ment; S is the shearing force; EI denotes rigidity (modulus of elasticity 
_ times moment of inertia); and p is the load per unit length on the beam. Both — 
i equations show that if the beam is loaded with the M/E I load instead of p the — ; 
differential relation between y, @, and M/E1I is similar to the differential 


relation between M, S, and ing tothe detintion 
_ This fact is the basis for the c conjugate beam according to the definition 


given above. It has to be remembered that ‘the conjugate beam is always 
statically determinate and each type of beam is associated with one whose 


_ “ELECTRICAL L BEAM” AND ITS ITS CONJUGATE E BBEAM 


 *Fig. -1(c) shows the complete circuit for the simulation of the simply sup- 
pected beam indicated Fig. 

3, each one of 

3 analogous to a length of the total lof the 


re “An Electric Model of Beams and Plates,” by E. Truds6, Byningstatiske Meddel- —_ 
elser, Vol. 26, No.1, August, 1955,p.1. 
“Theory of Modern Steel Structures,” ” by Grinter, r,Vol. T , The MacMillan Co. 
New York, 1949, p. 49, 
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fixed end in the principal beam becomes a free end in its conjugate 
beamandviceversa, 
ed by an un- 
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means of voltages ¥, measured with respect to points ‘Aor Bon the elec- 
tric beam, with voltmeter. Vj connected to either one of these points, and one 
of the accessible points (oft the string of of resistors. From these latter points, 
uniformly distributed on the beam, currents can be tapped by means of ree 
_ sistors R’, If these resistors are proportional tothe E I value, at each section 


= constant of proportionality). 


currents in the resistors F R’ is and assuming negligible 
voltages in the electrical conjugate 


A onstants c, and» ‘oF. related by! c’=—., Therefore, the 
currents in the circuit of Fig. 16) give values proportional to 
4 
To avoid appreciable change ‘in distributions (pro- 
portional to shears)2 in the electrical beam AB, currents I’ derived from the - 
principal beam must be small enough. This is obtained with values R’’ >> Ar 
a large value for constant c of Eq. 3. If a second electric beam P 
of the same Boats wat as the principal one is loaded with the currents I’, the - 
analog of the conjugate beam is obtained. = 
pes Currents in the conjugate beam will be proportional to the shears in it, 7 


which in turn are proportional _to slope in the elastic curve of the actual 


moe 
measurement of potentials along the beam will be done by. 
means ‘of voltmeter Vo with one of its terminals connected to a point of zero 
potential on the conjugate beam; i.e. to a simple support, a free endora — 
is connected to the simple support 
EXAMPLES OF TYPICAL INDETERMINATE BEAMS av 


In order to obtain the results in . the case of a determinate beam, such as . 
the beam in ‘Fig. itis enough to adjust resistor r Rp fixing 


ti ss The term Ap is an ammeter to measure the current Ip analogous to force © { ee 
ig ?P, which is fixed with the variable resistor Rp directly connected at Ap. — ae 
sa -_ one force is used as the load for reasons of simplicity. If more forces are ie 
- =_— they should be introduced into the analog in the same way as Ip. Cur-— = 
7 rents I, and Ip proportional to reactions Py and Pp are measured with 
— 

— 
— 
— 
— 
— 
® 


“then fixed. No more operations are > necessary. 


In the case of indeterminate beams the elastic properties of the beam must — 
also be considered. This can be done > by means of the conjugate beam, as will | 

_ be seen infourtypicalexamples. = 

The End Supported Cantilever.—Fig. 2 shows a supported cantilever 

beam conjugate beam (b), loaded with the moment diagram of the 

principal beam, and the electric circuit (c) for the simulation of both beams. © 

_Each string of resistors was replaced by one tapped resistor. This change can 

be made if enough taps are available. The power supply and the instruments 

- for measurements in the principal beam are not indicated since they are the © 

‘The conjugate beam, according to what was saidpreviously, has a free a 

at A’ anda simple support at B’. This is obtained in the circuit by leaving end 

- 3 ’ free and connecting B’ toB, the point of zero moment in the principal beam 

and therefore the point of reference for its voltages. 

a ‘The point ¢ of reference of the electrical conjugate beam | may be A’ because 
this is a free end and therefore it has zero moment in the mechanical beam, 

At and the point of A’B’ under consideration, 
oe _ Since the points A and B of supports have no caiman: the corresponding _ 

a A’ and B’ of the conjugate beam have zero moment. Voltmeter V2 when 

- connected between these two points will indicate zero. This | condition will | 
solve the indeterminancy of first grade for the beam. 

_ According to what was said above, the way of operation on of the circuitis 
follows: ~The c current Ip is adjusted to have the value to load 

P, then Ip is adjusted ‘until the voltmeter V2, connected to A’ and B’ ‘ indicates 
sero volts. Under these conditions currents and voltages on the electric beam 
AB will give -ne reactions and bending moments; voltages on electric beam 
AB’ will give the deflections of the supported cantilever. Iterations will be ° 
necessary if changes when Ipisadjusted. 

(b) The Fixed-Ended Beam.—The two indeterminacies existing inthis case 
are solved with the circuit of Fig. 3. Tare 

a: The conjugate beam must satisfy the following conditions: ‘first, the e ends» 
are free as indicated in Fig. 3(b), because in the principal beam they are 

fixed (note section entitled “The Conjugate Beam”); second, both ends have 


‘moment [elastic c curve of Fig. 3(a)]; third, the angle between n tangents | to 
the static curve at both ends is zero[Fig. 3(a)), 
a The first condition is satisfied if no connections are made at the ends of 
the electric ‘conjugate beam, ‘Nevertheless in Fig. 3(c a connection has been 

made from end B’ (or A’) to a point Z, of zero moment on the beam AB, This 
pau is necessary to give a return path for currents I’ to beam AB, but 
it does not constitute a “support” for the electric conjugate | beam because the — 


current at B’ has to. be equal to zero, which means zero reaction force and 
‘The second condition is obtained when the voltage between A’ and B’ is equal - 
. It has to be remembered, _ for the third condition, that the angle be- 

- tween tangents at two points of the elastic curve is equal to the area under _ 
a the M/E I eure between the two points (firs st moment area 1 theorem). 2 The | 


ee @ 3 the value of “force” Ip (afewiterations may be necessary if several resistors 7 aN 
: 
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a CONJUGATE BEAM 

and Z, is gusratienet to the areaunderthe M/E I diagram and it will be — 

when the angle between tangents at both ends is zero, = 

To adjust the circuit of Fig. 3(c) for the fixed-ended beam we begin by 7 

adjusting current lp according to ‘the given load. Next any point Z 
7 as zero potential) is connected to the ammeter and than I, and Ip are adjusted — 
until the voltmeter indicates zero. If the current Ig indicated by the ammeter — _ 

is zero the fixed-ended beam is correctly set up and a of — 


@ CONTINUOUS BEAM WITH THREE EE SPARS 


| 


FIG. 5. —CONTINUOUS BEAM — THREE SPANS, 


zero another Z potnt 
chosen and the procedure until voltmeter andammeter (two conditions 
for the two grades of give zero 


inacy ime zero moment at pen ‘end and, ‘since there a are no > deflections at the 
three go stand the corresponding ge A’, B’, and C’ of the conjugate beam 


= 4 
— 
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(c) ELECTRIC CIRCUIT FOR THE SIMULATION 
— 
a 


the circuit of Fig. ‘in n which the electric beam has been divided in 1 
parts, AB and BC, for each span. Point Z whose potential will be taken as 
reference, is directly connected to A and Cc so that th these points | will have zero 


= ends of the electric conjugate beam, A’ and C’, are coin ——€ 
and they are directly connected inthe circuit to the points / A and C of potential 
reference in the principal beam. Point Bt the conjugate beam is an unsup- 
ported hinge. In the electric circuit, B’ must satisfy the condition of zero , 
moment which can be obtained adjusting the reactions of the principal beam _ 
with resistor Rac until the voltmeter connected between A’ and B’ indicates 7 
Zero. Thus, when this condition for the indeterminacy is satisfied and currents | 
I, and Ig, representing the given forces P and P9, are adjusted, the circuit 
= is simulating the continuous beam of two spans; reactions are measured as 
the currents at points A, B and C; moments and deflections are measured as 
voltages on beams ABC and A’B’C’, using point Z as the zero potential ref- | 
— (@ Continuous Beam o of Three and More. Spans. _—Inthe case of three spans, 
¢ grades of indeterminacy have tobe considered. Fig. 5(c) gives the electric 
circuit for the beam of Fig. 5(a) and its conjugate, Fig. 5(b), 
rs, As in the previous case, ‘both ends A and D, having zero /moment and zero 0 
voltage are directly connected to each other and any one may be taken as the 7 
Zero potential reference, Points A’ and D’ are again connected to the ends A 
* and D and, since there are hinges at joints” B’ and C’ of the conjugate beam, 
these points must have zero voltage in the electric circuit; two voltmeters, 
connected between A’ and B’ and between B’ and C’, or D’ and C’, when indi- - 
cating zero voltage, give the two conditions re required | to solve — marcez: de- 


4 oetians must be varied until the voltmeters ‘indicate zero values. ae 
For each new span which is added, another indeterminacy is mbenteent x and 
_ therefore another simultaneous condition, for zero voltage between hinges of 
the new interior span, has to be satisfied. The process of iteration for ob- 
taining the simultaneous voltages on the conjugate elec electric beam 


‘cumbersome for. each’ indeterminacy 


‘SCALE 
: Moments are related to forces and seidiiiitaat by the equation M = F land 
__ voltages, according to Ohm’s Law, are related to currents and resistances by 
the e equation V = a. R. With these . relations between the 1e magnitudes ¢ of both - 
= (indicated in the section entitled “The Conjugate Beam”) the scale am 
factors for the principal beam can be fixed. It should be noticed that only two 
scales can be arbitrarily « chosen. ‘When this is done the third scale becomes © 
example, if 1 kip equivalent to 10m A is chosen as is the scale for forces: 
ie 1 ft equivalent to 10.0 as the scale for lengths, then  - 
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Therefore the scale 1 ft-kip equivalent to 0.1 volt slate from the choice of 


theothertwoscales, 
- The scale of length for the : conjugate beam may be the same as that used 
forthe principal beam, 

. ¢ In order to fix the scale of deflections that are measured by the voltages 
on the conjugate beam, the following data are necessary: the ae 
scale of the principal beam, the relationship between the E I values and re- 

= sistances R’ that join both beams, andthe Ar value corresponding to each one 
of the parts into which the beams w were divided (Fig. 1). eer ne hina 

_e An example will show how the scale for deflections is obtained. cas» 
- Consider a beam of 100 ft length with E I= 1010 tb in2, The beam is ai 
vided into ten parts, each 10 ft long, in order to load the conjugate beam with 
moment diagram of each one of these parts. The following scale relation- 


have been chosen. 


2. 1 ft ie : / 1 Volt (in the principal beam) =. 
3. 1010 Ib in? : 10°.2 = 
wih 

og If the beam were not of uniform cross-section or material the R resist- 
ances in each one of the ten sections of the beam would change. Each value of 

_ R’ would be fixed according to the E I value at each point as given by the last — 
_ Suppose a moment M= 1 ft kip of the p1 principal beam, along a 10 ft section, 
Ae loads the conjugate beam with its resultant M/E I diagram. In the analog a 


current” through R’ will the electrical conjugate sot that the fol- 


ee x x 10- 6 equivalent to I’ 

_ and since 1 ft is equivalent sali 


A 2 in 1. equivalent to19 


Multiplyin Eqs. 
Mu ultip ving Eqs q 
12,8 x 10°F in, 
Itages of the order of millivolts will be e measured in the 


electrical conjugate beam the final relationship used is: Tie 


(172, . x 10-3 in. deflection 1 equivalent to 1 mV 


| Sepertinedts we were made with an electrical beam simulating a a continuous 
beam of three equal Spans | of 100 ft each. A flexural rigidity EI= 1010 Ib in? 
_ The electrical beam was constructed with the resistance of 10002 for ea each : 


to represent a force of 1 kip. 
With these values the scale factors are those ‘computed in ‘the previous 
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CONJUGATE BEAM 


for the reaction at one end, the influence line for the moment at an inner — 
3 support and the influence line for the deflection at the enn of the outer 


, The experimental values can be directly compared to the corresponding 
given by the theoretically. calculated influence for the reaction 
_ Since the theoretical calculation of the entire influence line for deflections — 
would be extremely lengthy, only one point of this influence line has been 
computed: | the deflection of the midpoint of the first span when the load is i 
the same point. The resulting deflection was 2.52 in, = 
In the electrical conjugate beam under the analogous conditions a voltage 
of 14.5 mV was measured. This value multiplied by the factor 172.8 x 10-3 
in per mV gives a deflection of 2.50 in. Enough other deflections were ex- 
-perimentally determined s so influence line « can be readily dr drawn. = 


The. analog circuit for the simulation of bending of determinate 
and indeterminate beams by applying the method of the conjugate beam —_ : 
"mits the determination of magnitudes, including deflections, within the ac- | 
curacy of the electrical instruments and resistors used. The ease with which 
the electrical measurements may be made and with which the electrical 

- circuit can be changed makes this method of measurement very convenient — 
cases where immediate results are required. 
These electrical beams may also provide abasis 
the influence lines on Aa cathode oscilloscope or on an x-y recorder. 
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= constants of proportionality; 
= ‘flexural rigidity, ‘pounds x inchs 
= electric | currents, amperes; — 
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APPENDIX.—NOTATIONS 


P, Pi, Po, 


= load ‘per unit length, Kips per foot; 
Ra, electric resistances, ohms; 


wh 


al = = voltage on on the cir vo 


= 
ctrical conjugate beam circuit, volts; _ 
ion, feet, inch; and bac 
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_ METHOD OF INTERNAL CONSTRAINTS AND ITS APPLICATION a 


By Enrico Volterra,! 1 F. ASCE 


An that is of both and ‘systematic character 
presented. The method unifies a number of engineering theories in the statics _ 
and dynamics of plates and bars, and, in addition, is simple and direct. | —— 
method is based on assumptions concerning the components of the elastic dis- 
placement. In fact, it is assumed that the components of the elastic displace- 
ment must comply with special equations of constraint. - These equations of 
constraint, which determine the elastic displacement, are the basis for the — 
term “method of internal contraints. 
_ The method is applied first to the discussion of vibrations of elongated 
straight rods. The equations of motion and the boundary conditions for | flexural. 
vibrations are derived. _ It is shown that these equations reduce, in the plane- 
_ stress case, to the well-known Timoshenko beam equation. The equations of 
motion and boundary conditions for longitudional vibrations of elongated rods 
are derived, and the results given by these equations are compared with ~— 
_- by use of other engineering theories and with those given by the exact 
- Pochhammer solution. - Finally the method is applied to the discussion of flex-_ 
= ural vibrations of rectangular thin plates, and the equations of motion and pro- 
per boundary conditions are derived. These equations of motion coincide with 
the equations of motion derived by R. D. Mindlin, M. ASCE. _ The present equa- 


tions reduce, in the static case, to equations that are very similar to those ob- 


_ Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 


Society of Civil Engineers, Vol. 87, No. EM 4, August, 1961. 
- Prof. of Engrg. Mechanics, Univ. of Texas, Austin 12, Tex. 
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The classical theory of transverse vibrations of clastic. rectilinear bars, 
by Bernouilli and Euler, neglects the effects of rotatory inertia and of shear. 
If the results given by use of this elementary theory are compared with those 
_ given by use of the exact theory of L. Pochhammer and Chree (1)2 in the case a 
- of cylindrical bars of infinite length, it is seen that these two theories produce > 
considerably different values. An approximation due to Lord Rayleigh (2), : 
which takes into account the effect of rotatory inertia, does not improve the 
— -—In 1921 ‘Ss. Timoshenko (3) showed the important role played by shear in 
transverse vibrations of bars. | By adding the effect of shear to the effect of — 
rotatory inertia, he obtained a good approximation between the results given 


_byuseofthe exacttheory, 
Results as described inthe first paragraph of this Introduction are obtained _ 


j also in the case of longitudinal vibrations of cylindrical bars. _ Paras 
In the case of flexural vibrations of thin plates, too, there is no agreement 7 — 


exact theories of Lord Rayleigh (4) and H. Lamb (5). 

_ For static deflections of thin elastic plates G. Kirchhoff’s classical theory 
(6) neglects, in the study of deformations, the effects of shear and of normal 7 
stress. In many cases use of the Kirchhoff theory reaches results that are 
contradicted by experiments. For example, the application of the classical - 

_ theory to the case of a rectangular plate simply supported and submitted to a 
uniformly. distributed load on the free surface leads tothe false . result that four 
concentrated reactions must be applied on the four corners of the plate in order _ 

_ E. Reissner (8) (9) (10), H. Hencky (11), L. Bolle (12), M. Schafer (13), 

and A. Kromm (14) (15) have presented new theories of deflection of thin plates 4 
that take into account the | effect of shear. In these theories, whichare referred — 

to in the text as engineering theories in order to distinguish them from the : 

exact theories obtained by applying the equations of the mathematical theory 
(Of elasticity, one generally starts from hypotheses" concerning tl the stress dis- 

tribution in the plate in order to reach the equations of motion or the equations 

‘For example, in his theory Reissner starts } from appropriate and rational 
hypotheses concerning the stress distribution across the thickness of the plate. 7 

He supposes that the stresses due tothe bending moment vary linearly, whereas 
those due to shear vary according to a parabolic law across the section of the - _q 
plate. At the same time Reissner does not neglect the effects due to the normal > 

Stress. Reissner has shown that his of equilibrium satisfy three 
boundary conditions, whereas” the equations established by Kirchhoff satisfy 
only two boundary conditions. . The equations | of equilibrium obtained (8) by 


a Reissner in 1944, taking into account the transverse deformation of the plate 


RET ESO 
— 
tained by EF PReiccner Unwever the nrecent method of deri the equations 
is simpl 
x 
| 
ae 
a 
ae 
- 
Ps 2 Numerals in parentheses refer to corresponding items in Appendix II. “eq? ae = 


to shear by applying Castigliano’s theorem, were derived directly by A. E. 
7 By applying his S equations Reissner was able to solve various problems of 7 
practical importance, including the problem of bending and torsion of an infinite | 7 
a ‘Plate with a circular hole in it, and the problem of bending ¢ of a cantilever plate 7 
’ _ due to transverse forces applied at its extremity. In the foregoing cases the 
oe application of the Kirchhoff classical theory leads to results in contradiction © 
with» ‘reality, whereas the application of Reissner’s equations of equilibrium | 
leads to results which are in accord with experiments. 
_ Bolle, independently from Reissner, obtained the same equilibrium equa-— 
vs tions and boundary conditions in his important work on the deformation of thin 
plates taking into account the effect of shear (12). ee eS =. 
oe In this paper the writer presents a method that is a approximate and at the 
same time has a generaland systematic character. The method unifies a cer- 
tain number of engineering theories in the statics and dynamics of plates and — 
bars. It is, in general,a more simple and direct method. Whereas the methods 
_ mentioned previously are based on hypotheses ‘concerning the distribution of 
stresses, the method presented herein is 8 based on assumptions concerning the 


 ponests of the elastic displacement must ‘comply with special equations a 
ae These conditions of constraint, which determine the elastic dis- 
placement, are the basis for term “method internal constraints. 


_ plane-stress cane, to the well-known Timoshenko beam equation. The equa- - 
tions of motion and the boundary conditions for lc longitudinal vibrations of elon- 
i "gated rods are then derived, and the results given by these equations are com- 
_ Pared with those given by use of other engineering theories and with those given 
i“? use of other engineering theories and with those given by the exact Poch- 
7 hammer solution. Finally the method is applied to the discussion of flexural 
_ vibrations of rectangular thin plates, and the equations of motion and — 


boundary conditions are derived. These equations of motion coincide with the 


equations of motion derived by Mindlin (17). These present equations Teduce, 


vibrations are then derived. It is shown that these equations reduce, the 


= and more = 


_ FLEXURAL VIBRATI ONS OF STRAIGHT ELASTIC BARS. 


_ Netation.— —The —The letter symbols used in this paper are defined where they first 
appear in the text, and are arranged alphabetically for convenience of — 


ence inAppendixl. 


For flexural vibrations the components of the elastic displacement u, and 
w in the direction of the axes are assumed tobe: 


X, Y, Z; t)=ya t 
u( 


9 y2 
t) = vy teu (x; 


| 
| — 
1 
> 
| 
ee 
a 
— 
present method of deriving the equations is sim- 
Paley 
— 
(Ib 


 t)=y2fy t) 


(Concerning the geometric of the elastic constraints 
by Eqs. 1, 14, and 24 reference is made to a previous paper (20)). In Eqs. 1 


- is time and is Lame’s constant. The axis X coincides with the axis of the 
‘bar, whereas the Y and Z axes coincide with the principal axes of yet of wall 


Oy = ( +! 


026860 


Eqs. 2 show that, in order t to satisfy tl the conditions that = = 


lowing relationships must be verified: 


L 

~ 5 al 6, terme. containing y2 and y z have been neglected 
inthe expressions forthe components of strain and particle velocities in com- | 


parison with linear terms. The term L L Tepresents the length of the bar and 


_tangular of breadth Band height q 
m= j. From Eqs. 1 the following expressions are found for the normal stresses _ 4 a As 
ma 
which and G (the modulus of — 
O,the fol- 
_ From Eqs, potential 
fi 


INTERNAL CONST! CONSTRAINTS 


‘Se the moments of inertia of its cross-section with respect tothe y y and z axes, 
and A the area of the cross-section. The factor k is introduced into Eq. 5 to 
take into account the non-linear distribution of shear- stresses | over the cross-— 
‘section of the bar. It has the same significance as that of the well-known Tim- 


“oshenko shear coefficient. 


On the that no external rinciple t (body | and surface forces) 


should be ‘pointed out that the validity of 


cr and 8 is not limited to “a 
In fact the only parameter related to the 
tion is the ratio of the moments of inertia ly/lz, so that 
7 and 8 2 are fora bar of any cross-section, 


By eliminating th the function ax in 7, they can “i under the fol fol- 


lowing form: 
2 tee 


“bar of cross-section. 


+v) 


—~ 
— 
im 
a 
— 
2 
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i In Ba. 11, 2 the terms due to shear are neglected, the differential equation - 
ist modified to the form usually attributed to Rayleigh: aiid 


por to the following « equation: 


of the elastic are ‘considered, the following equations mo- 


tion is | obtained: 


2 


"From Eq. 10 the equation is 


“4 2 


the bar 


4 
es 
— 
= 
wa 
} and to shear are neglected, the 
— = im 
iy which A is a eonstant i= -~-1 ~= — “Jf, is the wave-lencth and c the phase 
v 
= 
= 


INTERNAL L CONSTRAINTS | 


In Fi g. 1 phase of flexural waves in cylindrical ‘bars according 
to the elementary Caeory, the Rayleigh the theory, the Timoshenko theory, the ex- z 


Elementary 


s 

Method all Interna 


al Co nstraints 
ion) 


(First “Approxima ation 


Timoshenko's Theory 9) 


4 Diameter of the bor 


_ Wave- e-length 


BE 


PHASE VELOCITY OF FLEXURAL WAVES 
IN CYLINDRICAL — 
act Pochhammer theory, oh and the present method (first and second ap- 
are compared. 


the case of in a bar cross- -section 
the components of the elastic displacement will be assumed to be: oo 


X, y, Z; t) = vy (x; 


4 
a 
| 0.5 Exact Theory according to Pochhemmer 
“ae 


, 1961 


23 (x; £3 (x; t) {4 (x; t t) 
ae. The om a f1, fo, f3 an and f4 are e determined in terms of the einen: func 
tions Vx; » Ay and Az from the condition that the free ‘surfaces of the bar are 
free from normal stresses. 


= 0 tor y 


+ 


variables y ‘aad z, they are equivalent to the following four equations: | 


4 2 + (+ + (A + dy 


(142.6) 


9g 


y 


| 
H2 


In the evaluation of the potential and kinetic energies the er, terms con- 
taining r 22, yt and z4 are neglected in comparison with second degree terms. 


;t)..(17>) 
|: 

m (A+2G) \ay+- z2f3+—-f4)=0.. (18) | 
| 
18 and 19 need to he satisfied identically for anv value of the 
(20a) 
— 
+ (1 +2.G) ay | 


INTERNAL CONSTRAINTS > 
ct to » variables y and z the following eupressions 


\ 

| 


In the absence of body and surface forces, applying the Hamilton on principle 


4 


a The equations of motion and corresponding boundary cc conditions for the | plane : 
stress ¢ case are from Eqs. 23 and 24, ‘In fa by letting B approach 


Lit 


For the at =0,x=1: 


Ay | 


— 


length c the velocity, from Eqs. 23 the equation for 


> 


V 


te 


24 - 
The constant t Bi in Eq. (28d ratio the sides ‘B/E. 


jj Substituting this value for Into Eqs. 23 and 2% the equations of motion 
_ and the boundary conditions in the case of plane stress are obtained — relies 
(a) For the equations of motion: | Wa * 
q 

= 
— 
— 
> 


-0 in Eq. 28d the the case of of plane stress is 


case of a bar of square cross-section = 
with the results obtained by the use of other satiaoamtee theories of longitudinal 
wave propagation and with the results obtained for a oes bar | from the 


cular cross section and square cross section because, according to 
made by R. W. Morse, M. ASCE (22)(23), the dispersion curve for a bar of 7 

_ square cross section is very close. to the one obtained for a bar of circular 

cross section if the ratio of the diameter to the side of the square cross sec- 


‘The engineering theories ( results of use of which are compared in Fig.2 ff 
with results obtained by use of the second approximation of the method of in- © 
ternal constraints) are the elementary theory (curve 2) (24), Love’s theory 
(curve 6) (25), Bishop’s theory (curve 3), Bishop’s improved theory (curve 4), 
_ Love’ s theory improved by Bishop (curve 7)(26) and the method of internal - 
_ constraints (first approximation) (curve 5). The curves for the exact theory 
according to Pochhammer for cylindrical bars (curve 1) and the method of in- 
ternal constraints (second (curve 8) are also shown. Because 


derived, and their implications “discussed, by the writer in a previous paper 
(27) they will not be discussed here. _ 


VIBRATIONS OF THIN PLATES 
Inthe case of flexural vibrations of a thin —— plate the components 


the elastic displacement will be assumed to be: 


= Ux (x, y; t t) (30a) 

t = Z My (x, - (30b) 


- The rectangular iia of sides a and b is referred to a system of orthogonal 
- cartesian coordinates X, Y, Z, the plane X Y coinciding with the middle surface 
of the plate and the Z axis | directed downward. 4 


— 
— 
— 
q 
&g ae 
4 
— 
.. 


~ 


r radius of gyration of the cross-section of the bor __ 
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FIG, 2. —PHASE VELOCITY. OF ‘LONGITUDINAL wes 
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“components of strain computed fr from n Eqs. 30 are: 


OZ My + ay ay 


= 


tained from Eq. 32 for the unknown yit t): 


>= 


“order terms z2 will be Ye neglected in ‘comparison with linear terms. “After. 


4-4 tegration with respect tothe variable z the following expressions =e obtained: — 


(a) the energy: fi... 


Ux \2 vE = 8 Ux x 2 by 
a 2 1-» 


oy 


Oz = 1+ 2G eg = Az poy * +2Gzf ...... (32). 
— 
+ 2 k H Ux Ox | \ dx dy ( 34 


‘Fort the kinetic ic energy: 


* 
in which Hi is the constant ‘iia of the plate, and A is the area of the sur- — q 


face of the e plate. factor k (shear coefficient) was introduced, in the 


previous cases of the bar, to take into account the non-linear distribution of a 
% the shear over the cross-section of the plate. 2s 


names Hereafter the function Tu, which ‘represents the vertical displacements of 
a middle surface of the plate, will be indicated by the letter Wo. 


The work by the external distributed for force p(x, t) acting on the 


‘By applying Hamilton’ s s principle, the following equations 0 of motion: 


= 


(38a). 


4 


nas 


— 


we = px, y; 37) 
, A; 
— 
4 
— 


‘In Eqs. 38 and 3 nd 39 


The shear forces Tx and Ty, the bending moments Mx and M ee twist-— 
moments S Mxy and Myx per unit of the defined by 
the following expressions: 


, Oy z1dz= + 


» 
Ty * — 
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Invi w of 42, and 44 4, Eqs. 38 can be w on: 


May 3 Pig 

x ~ 

_ 


3 The t terms on the right of Eqs. 45 represent the terms due to rotatory in 


ertia. Eqs. 45 arethe same equations as those obtained in 1950 ina more 2 
plicated way by Mindlin (17). 


il Differentiating Eq. 38a with x respect to x and Eq. 38b with respect to y, and ~ 


. 


if the terms due to shear a are neglected in ‘Eq. 47 it becomes 
4 


(os 
D4- 

ne 


4 y,t) + ‘pH 
if terms due rotatory inertia are neglected i in Eq. 47 it 


= p(x, yt) .. 
Vom 


4 p (x, y, t) 4 


lected at the same time, it reduces to the classic plate equation 


im 
(45>) 
4 
eliminating the term + —HY between Eq. 46 
— 
4s) 


wae In Eqs. 2 components of the elastic displacement, if only the first 


terms are considered, or if: 


u(x, y, 2, t) = = (x, y,t 


7 


v(x, Z, ,t) =2 = 
xX, 


hen the eae equation of motion is obtained in final ssakatl 


x, y, t) 1- 


in which A is a constant, i = ~— 2 1 /L, L is the wave length and c the 
_ From the equations of motion (Eqs. 47, 48, , 49, 50 and 52) the following dis 


(a) From at 


c 
= 


For 
Co/ | 


EM4 
D Ad wo (x,y, t) + pH = p(x, y,t) .....(50) 
— 
4 
x,y,t) 
7 
(= a 
q q 
| From Eq. 48 


co) 
In Fig. 4 the curves c/co) against H/L, according to es: “54, 55, 56 and — 


31, are suageens with the results obtained by use of Lamb’s exact theory . The 


|.Elementary theory 2. Method of |.C. neglecting shear 

_ 3. Method of 1.C. Ist approximation 4. Method of |.C. 

neglecting rotatory inertia 5.Exact theory according 

ofl. C. including shear and rotatory inertia 


= 


4 St 
| 
| 
| — 
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DEFLECTIONS OF THIN RECTANGULAR PLATES 


In the static case 


- 6x 


7 while the boundary conditions remain th the same Eqs. 


+ 


The equations obtained from 42, 43, 44 60. 
- 


q 
— 
by 2 ~ 4.2 2(1+v) by 2(1+v 
while Eq. 46 
7 
...(62a) 


“In the case a very thin and 61, 62, 63 reduce 
the equations ate the classical 1 theory: 


pa A a Wo=P x, y (65) 


| 


D— 

‘Eqs. 61, 6 62, 63, and 64 conreapond to the following equations of Reissner’s 


Pp 
10 = 1-v) ‘x2 


a. 
+ 


thee case the value of the shear coefficient k should be k = 1/_ 

q 

— 

> 

p 66 ( 69) 


Mxy 10 \ dy 


2 can shown that, theory th that is presented, as with the 
a theories of Bolle and Reissner, the solution of equilibrium problems for thin. 


(73) 


that defines the stress function 
Shear forces are computed from the 


iT 


CONCLUSIONS 

— The approximate method sennsalad is applied to studying, in a systematic 

ye! and general way, static and dynamic problems of plates and bars, and it unifies 


of engineering theories on the subject. = 


By applying the method of internal constraints (which is based on special 


assumptions concerning the components of the elastic displacement) , the equa- 
- tions of motion or the equations of equilibrium and the corresponding boundary 


conditions are derived more directly and more simply than through application 


f engineering methods based, in general, on hypotheses concerning the a * 
bution of stresses. 


‘“ In the case of flexural vibrations of elongated straight rods, the ‘equations — 
fe) 


f motion are e reduced in the p plane- -stress. case to the well-known Timoshenko > 
aa Inthe case of ' longitudinal vibrations | of straight 1 rods, when dispersion curves — 


ed 
4 + Ky 
| 
| 
— 


obtained by use of other engineering methods and by use of the exact Pochham- a 
mer solution, it is seenthat the former curves give results closer to those | ob- : 


tained from the exact solution. 


- In the case of flexural vibrations of rectangular thin plates, the equatio 
a motion derived by applying the method of internal constraints coincide with 

ae = equations of motion derived by Mindlin from the three-dimensional equa- _ 
tions of elasticity. However, the present equations are obtained in a more di- | 
- way. These - equations are reduced, , in the static case, to equations that 
are very similar to those obtained by Reissner in his theory of bending of thin 
taking into account the transverse-shear deformation of the plates. 
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APPENDIX L- “NOTATIONS 


The following symbols are used in the text: 


= _cross- -section area of beam; 


= breadth of the beam; 


= velocity of propagation of waves in general; 
=~ 


velocity of longitudinal wave waves s of infinite -lengthin 


4, ‘fy, fo, f3, which | define the elastic 


modulus of rity, G: = 


of {inertia of a area about the Y and Z 


= twisting moments; 
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external distribute 


= radius of gyration; 


= time; 
= shear shear forces; 


= kinetic energy; 


« in general; 


rectangular axes of ne 


= rectangular coordinates; 


= shearing strain in rectangular coordinates; | 
= unit elongations in the and Y, and Z directions; 


> 


Txys Tyzs shearing- stress components in rectangular’ coordi- 
nates; and 
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7 aoe 


COMPARATIVE STU STUDY OF A SEGMENTAL ARCH RING! 


C. ZIENKIEWICZ ASCE. —It is gratifying th that the co reached 
by A.  Veltrop, ASCE, with the wide range | of values investigated ve ify 


7 _ While the object of the paper was a comparison of various shell theories on 


elementary example, it is interesting to see some of the results including 
the elasticity of the abutments placed side by side w with these. It is rather re-— 
-markable that the effect of this is of the same order as the discrepancies be- 

Bee the results of the various theories applied | to the rigid abutment exam- = 
ple and points to the caution necessary in interpretation of such results, - 


| 


January, 1960, by O. C. Zienkiewicz (Proc. Paper 2332). 
Prof. of Civ. Engrg., Northwestern Univ., Evanston, 


— 
3 
— 

A 


— 
So 


Ve 


GEORGE HERRMANN NN,6A ANDAN} ANTHONY AR ARMENAKAS, 7 MEMBERS, - 
The writers concur with ‘Masur’ 8 view that space might have been conserved if 
the results were presented in a more compact notation. However, they do not 
share his opinion that the use of cartesian components resulted in a measure 
of obscurity, for if that were the case, this notation would not have been used a A 
by the many investigators cited in the bibliography. Masur is correct in stat- 
ing that the matrix 7;; does not transform as atensor evidenced from its defi- 
nition and fromthe tensorial properties of Ojj- . The writers, however, , disagree 
» with the suggestion that the principles of thermodynamics limit the choice of 
e acceptable stress-strain relations to those resulting from Eqs. — 
“\ The formulation of the problem considered involves the assumption that the 
influence on the motion of the components of strain and rotation, associated 
with initial stresses, is negligible. In a sense, therefore, the initial state of 
stress may be conceived as a state without any related strain, Hence, any ap- 

i proach which associates a strain energy density function with both the initial © 
and final state of stress may be misleading, 
3 In the theory presented by the writers, the additional elongations, shears | 
and rotations are smallas compared to unity, therefore, it is assumed that the 
difference between the initial and final stresses is related to the additional 
strains: in accordance to Hooke’s law for an isotropic | elastic solid. However, 

7 a” in the presence of initial stress (even in this range of small deformations), 

- the Trefftz components of stress 03; (whichinthe vocabulary of the mathemati- 
-_ cal theory of finite elasticity are referred to as material stress components) — r 
and the actual components of stress 74; (defined in the original paper) cannot 
be considered identical. Thus, the presence of a number of terms involving _ 

initial stresses in the displacement equations of motion is dependent on whether _ ) 


‘accordance to Hooke’ s law. It is this point which ‘the: writers attemptedto em-— 
 phasize inthe original paper by presenting three different possibilities for in- — 
_troducing the initial and additional stresses into the non-linear e equations of. 
T motion. Inother words, a number of terms involving products of initial stresses a 


June, 1960, by Herrmann in and Anthony Armenakas (Proc. Paper 2500). 
6 Assoc. Prof. of Civ. Engrg., Inst. of Flight Structures, Columbia Univ., New Touk, 
7 | Prof. of Civ. Engrg., The Cooper Unior for the of Science 
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order of magnitude as terms disc 
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= 


¢ the principle of adele the maximum resisting moment of a member, 

using the stress-strain relationship of Smith and Young are presented, | The 

same approach was proposedin 1957 by another writer18 who useda parabolic 


_ stress-strain relationship of concrete to demonstrate the principle, 


a 
The -_stress- -strain relationship, proposed by the authors, is unsuitable for 
design purposes because of its complexity, and it was developed primarily for 
= use in research, It was applied to the case of over-reinforced beams to show 
-thata simple ultimate strength design formula may be arrived at by analytical 
than experimental procedures. 


The e knowledge of the stress- -strain re-— 
~ lationship also makes it possible to predict the behavior of a reinforced con- 
crete member throughout the loading range with accuracy, as was recently 


-stress- strain relationship for various concretes at different rates: of loading ‘ 
field. 


_ The authors agree fully with Ali’ s suggestion that determination of the 7 
- ‘desirable. A p paper of Hubert Rtisch20 published i in 1960 reports « on work in 


e, 1960, ic Ladislaw B. Kriz and Seng-Lip Lee (Proc. Paper 2502), 
16 Assoc. Development Engr., Structural Development Sect., Portland Cement Asson. 


_ 17 Prof. of Civ. Engrg., Northwestern Univ., Evanston, Ill 

«*'18 “Het Gedrag Van Betonconstructies Onder Invloed Van Kortstondige Belastingen, ” 

IBC Medelingen, Vol. 5, No. 23rd, Delft, April, July, 1957, pp. 52-62, 
419 “Ultimate Strength of Non-Rectangular Structural Concrete Members,” by A. H 
: Mattock and L. B. Kriz, ACI Journal, Vol. 32, No. 7, , January, 1961 (Prossssings Vol 
57), pp. 737-766; Bulletin No. D48, PCA Development Dept. | 
20 “Researches Toward a General Flexural Theory for Structural Concrete,” 
_ Hubert Rusch, ACI Journal, Wel. 32, No. 1960 Vol. 57), 1), PP. 
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Closure by Morris Ojalvo 


is afurther development of the of other ¥ writers. His ‘work in 
"stems from Theodor von Karman, Hon. M. ASCE, who recognized that the beam- 
column was a stability rather than abuckling problem and Chwalla2 who appar- — 
- ently was the first to recognize that segments of the column deflection curves’ 
= were equivalent to te equilibrium configurations of a column under © 
increasing load or ‘moment. The outline for determining the limit of stable 
Tg ~ satay for any condition of end moments and restraints and the manner 
in which the information from the column deflection curves can be arranged 
> in nomographs to give rapid solutions are the only new contributions of of the 

we It was not the intention of the writer to imply that a compressed member 
occurring in a structure with sidesway could not be considered as a segment 
of a column deflection curve. This" aspect was not investigated because the | 


7 _ analysis of frames with sidesway was not the subject of the paper. 


= to rapid solution by nomographs. The principle involved, however, ae 
same. One determines successive equilibrium configurations of the struc- 
ture under increasing proportional load until an equilibrium configuration can © 
: no longer befound, The maximum load for which an equilibrium configuration | 
can be found is the collapse load, 
fant proof is adequate to show that a differential equation (Eq. 15) 
_ describes the shape of the deflected member. _ It would have been more exact 
-_ however to state that in Eq. 15 the abscissas are measured parallel to the x 
axis as it appears in Fig. 9, and that y is the distance between the thrust line 


and the curve measured parallel to the y axis. Inthis manner the usage of the 


term abscissa would also correspond to its definition as foundin a mathematics 
_ dictionary.16 The assumption that the shear must be small compared to the 
: axial force (the respective components of the total thrust normal and parallel — 
_ to the chord connecting the ends of a member) is nevertheless necessary if a 
compressed member is to be considered as part of a columndeflection curve. 
os _ The proof that Eq. 15 is a symmetric periodic wave under the assumed 
‘moment- curvature relationship is most welcome for it emphasizes an erron- 
-eous made by the in the text. ‘The writer’s statement (first 7 


October, 1960, by Morris Ojalvo (Proc. ‘Paper 2615). 
15 Assoc. Prof., Dept. of Civ. Engrg., ‘Ohio State Univ., Columbus, Ohio. 
“Mathematics Dictionary,” by Armen A, Alchian, G. James and R. 
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in entitled Theory) that the entire column deflection 
curve may be constructed from a known quarter wave length should have been 
4 qualified so that the statement was applicable on only to columns of symmetric 
 eross section and of a material having the same properties in both tension and 
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GONDIKAS, 26 and ‘MARIO G. G. ‘SALVADORI,27 F,  ASCE.- —The writers feel 


produced “discussion by Toebes, because their ‘interest in n dynamics ‘has 
_ a them a clear picture of the need for more abundant information inthe — 

field of dynamic inputs for civil engineering structures. 

Toebes is interested primarily indynamic inputs, whereas the engineer 
is mostly interested in dynamic outputs. — The applied mechanician who tries 


"simple unrealistic. inputs. ‘but of being lost when trying to deal with ‘complex 
ealistic inputs that he cannot easily feed into his equations, its 
ae More realistic static wind loads of the type used by Placido Cicala, 17 are 
being used to obtain membrane stress distributions in shells. _ Unfortunately, 
shells will withstand almost any load by means of membrane stresses but are 
liable to give in at particular spots due to bending st stresses, The writers have — 
therefore considered of some use the first solution to the bending problem of 

‘shells under conventional wind loads and notice that the: same problem has been 

successively solved by a number of other authors, _ 

The expansion of any ‘wind loads into a double in @ and 

been obtained and has made possible, through the work of Gondikas andothers, 

_ the bending analysis of shells under anti-symmetrical loads. _ It is, therefore, 
foreseeable that the evaluation of stresses due to fairly complicated static — 
will soon appear inthe literature, 

In order to be able to attack the corresponding dynamic problem, the dyna- 
mic wind loads must be known with much greater accuracy before the struc- 

_ tural engineer will bother to usethem in a dynamic stress analysis. There is 

no doubt that what has been done for earthquakes will eventually be done for _ 7 

hurricane winds, The delay in the study of dynamic w wind stresses is probably — 

due totheir lesser importance and tothe large amount of ignorance in connec-_ 
tion with vibrations of complex shell structures, 


e 4 October, 1960, by P. Gondikas and Mario G. Salvadori (Proc. ee 2616). 
26 Partner, Gondikas and Righidis, Athens, Greece. 
_ 27 Prof. of Civ. Engrg., Columbia Univ., New York, ‘N. Y., also Paul Weidlinger, 
Cons. New York, ¥. 
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% Discussion by | G. 2. Reynolds 
™ 
G Cc. REYNOLDS. author’s use of an n equivalent E to se non- 
linear beams appears to be an unnecessary complication. In linear beams the | 
curvature ¢ is equal to M/EI and the curvature, in the form M/EI, is used to © 
a slopes and deflections. For the analysis of non-linear beams, the au- 
thor has replaced E by an equivalent value E, suchthat o = M/EL to represent 
curvature, A similar device, which has been ‘used, 13 isto replace M by M such 
7 that ¢ = M/EI. As far asthis writer is aware, no one has replaced I by I, such | 
that = M/ET. There is no point in using an equivalent E, M, or I, because 
each case, it is simply the curvature which is used. It is a simple matter to 
obtain a moment-curvature relation for non-linear beams!4 without introduc- 
ing an and the analysis of the beams is notas difficult as is com- 
‘monly supposed.!4 The curvature, obtained directly from the author’s 
20), is nt thle which e is the strain at distance y from the neutral axis. The 


a 2 December, 1960, by Annabel L. Tong (Proc. Paper 2690), M, 
12 Senior Lecturer in Civ. Engrg., Univ. of Tasmania, Hobart, Tasmania. by aw 1a 

13 “The Elasto-Plastic Analysis of Two Portal Frames,” J. 


Limit,” by G. c. Reynolds, Civ. Engrg. Transactions, Inst. of of Engineers, Australia, 
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by | Shan S, Kuo M. 


8. KUO, 9,13 M. ASCE, — The flexibility influence coefficient matrix [F] 7 
of a shear building has been. presented from the given | geometry y of the building. 
The standard matrix iteration procedure is then used by the author to find the 
frequency and the corresponding normal mode 
_ Once the matrix [H] = [F] [M] is known, the basic problem is to find all nu 
= eigenvalues and eigenvectors of the matrix [H]. This canbe accomplished - 
by making use of Jacobi’s method. 14 This method reduces a symmetric ma- 
trix to diagonal form by performing a sequence of plane rotations, eachof which — 
reduces to zero one of the nondiagonal elements of the matrix, until the sum of 
Squares of nondiagonal elements fails to decrease. This method can be opel 


vantage of this method lies inthe fact that a complete listing of the poner nel 
and the corresponding eigenvectors is available. An analysis of relation be- : 
tween order “ matrix and number of plane rotation has previously been given 
by the writer.15 While Jacobi’s method was of theoretical interest only when. 


~ 


s was presented in 1846, it has become quite popular with the advent of high- - 

‘speed digital computers, — 


Ste. 


4 February, 1961, by M. L. Pei (Proc. Paper 2728), its” 
13 Assoc. Prof. of Civ. Engrg., Tufts Univ., Medford, Mass. =~ 
14 “Uber ein leichtes Verfahren, die in der Theorie der Sakularstorungen vorkom- 
‘menden Gleichungen numerisch aufsulosen, ” by C. G. J. _ Jacobi, Z.reine angew Math., 
Vol. 30, 1846, pp. 51-95, 


ae 15 “A Note on Jacobi’s Method for Real Symmetric Matrices,” by S. S. Kuo, Journal 


- of Aero-Space Sciences, Vol. 28, No. 3, March, 1961. 
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by ALA AL 


A, A, EREMIN 11 M. ASCE.—Precision of dynamic stress computation va vary 


with the assumptions adopted in expressing the loading and elastic properties — 
eh The assumptions adopted by the authors requires further clarification, Ob- — 
viously, the finite method computation has an important advantage, that is it | 
may conveniently be performed by means of a computing machine, _ However, 
the basic characteristic of the finite method, when applied to a structure, is 
the continuity of structure material functions. In the simple span beam con- 
sidered by the authors, plastic hinges may be formed at the plastic stresses 
which require special consideration when the finite method of computation is 
| 
_ In applying the finite method of computation to the transient nt loading, the au- 
considered the interval of loading designated by k. 


_ In the illustrative example the interval of transient loading was assumed | 


5 ‘equal to one half of natural period of — 


‘The natural period of vibration is yetersing to the vibration within 
_ the elastic stresses limit. At the plastic stresses, the periodof natural vibra- 

shall be considered with limitations. In addition, in the loading expressed ; 

; in Fig. 1, it is hardly necessary to refer to the period of natural vibeation, 


 @ February, 1961, by Melvin L. Baron, Hans H. Bleich, and Paul Weidlinger (Proc. 
Cons. Structural ‘Engr., Berkeley, Calif. 
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4 MOIRE FRINGES AS A MEANS OF F ANALYZING ; STRAINS@ 
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Discussion by Stanley Morse 


that the authors have not introduced any new method of measurement by which 
the partial derivatives of the displacements can be determined with greater 


precision than by n methods | previously reported, Both the “seometrical” ap- 


MORSE, 13_ order to confusion, it should be made clear 
14 jointly by the authors the writer, and the 


least, finite deformations. . It is a step in the right direction for the analytical 
_ approach to ‘consider the tensor strain equations, as the authors have done, | 
: However, the distinction between the Eulerian description of strain, which the 


: - terns ‘tends to be limited to what should probably be considered large, or at 


moire patterns yield, and the Lagrangian description, on which the usual engi- 
neering strain is based, should not be overlooked. Also, caution is indicated 
¢ the application of approximation procedures, in’ which it should be mentioned 7 
that the series expansion in Eq. 19 is incorrect. te This should read, including — 
only: terms up to 4th powers 


‘to 


all be difficult to maintain that the method is either easier or ‘more | accurate 
than the ‘ “geometrical” one presented in the earlier paper.14 
Ingenuity n may find other possibilities but it seems that, on lines througha 7 
- given point in a_ moire pattern, there are four conceivably significant meas- 
urements. These are (1) fringe spacing perpendicular to the grid lines, (2) 
fringe spacing parallel to the grid lines, (3) fringe spacing along the common : 
_ normal to the fringes, and (4) orientation of the fringes with respect to some 
fixed direction. From a pair of photographs of moire patternsof models simi-_ | 


_ 4 February, 1961, by C. A. Sciammarella and A. J. Durelli (Proc. Paper 2736). 
13 Assoc. Engr., Mechanics Research Div., Armour Research Penn, Chicago, 
14 “Geometry of Moire Fringes in Strain by S. Morse, A. J. Durelli, 
4 -C.A.Sciammarella, Proceedings, ASCE, Vol. 86, No. EM 4, August, 1960. it~” 

15 “Utilisation des Reseaux pour l’etude des Deformations, by Dantu, Laboratot re 

"Central des Ponts et Chaussees, Paris, 57- -6, ‘1957. 
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It will be noted that Eq. 34 contains a 3rd power imaicating addi- 
a 1al consideration should be givento the criteria by whichterms in the series aes 
may be considered negligible. 
The problems connected with performing a practicalstrain analysis by use 
= 


grid lines, there will be a total of eight such measurements. Of these, in the 
application o of the large deformation theory, t the methods of the present paper 


- graph and, except on lines of symmetry, both photographs are required. 7 
The “geometrical” method, on the other hand, permits, from one photograph 
_ the determination of the strain perpendicular tothe grid lines and the rotation, | 
both inthree independent manners. To any claim that, in some cases, certain 
of these measurements | ((1) and (2)) cannot be made with accuracy, it must — : 
be said that any such inaccuracy reflects . directly on on the methods of the p pre- 
_ sent paper and that the “geometrical” method will always offer at least one al- 
ternative method of measurement. 
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larly loaded and similar in all other respects except for the direction of the q a 


_wu Lu,9 A.M. ASCE,—The authors have not used the 1 method of three- 
: -moment equation to obtain solutions to the problems treated in the paper. The 
well- -known| three- -moment equation was firstd derived | by F. Bleich10 (it is called 


succeeded the basic equation to include the effect 
of lateral loads applied between the joints. By utilizing the results contained 
in the present paper, it is possible to further generalize the equation so ir 
the effect of transverse foundation reactions can be taken into consideration, 
_ The three-moment equation in its most general form will be derived from 
the slope-deflection expressions of Eqs. 18(a) and 18(b). In Fig. 4 the de- 
formed shape of two consecutive members Lj and Lj+j of a continuous beam- 
- column acted on by the compressive forces Nj and Nj+1 and the transverse 
loads pj(x) and pj+1(x) is illustrated, The deflections yj_4, yj and yj+1, and 


~ @ April, 1961, by S. L. Lee, T. M. Wang, and J. S. Kao (Proc. Paper 2801). —- 
_ 9 Research Assoc., Fritz Engrg. Lab., Lehigh Univ., Bethlehem, Pa. | 
10 “Die Knickfesigkeit elastischer Stabverbindungen, ” by F. Bleich, Der — 
10, 1919, pp. 27, 71, 117, and 163. 
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CONTINUOUS BEAM-COLUMNS ON ELASTIC FOUNDATION@A 
the stability of plane frameworks. In its original form, the equation is only 
f 
7 
if 


the end moments Mj- My and will be taken as positive when they 

_ inthe directions as shown inthefigure. By applying the slope- deflection equa- 

tions tothe span Lj andeliminating the slope 6 j-1 fromthe two resulting 


tions, the following expression for the slope may be 


- 


cos hy sin ¢ - ¢sinhycosd 


g” + Sin hw o+ cos hy sin ne) 


ver 
Repeat the same process for ‘the span Li+1 and eliminate the slope Oi+15 the 


| The condition of continuity at the center a requires that the value of a 
_ given by Eq. 57 should be equal to that given by —_—s 59, By equating Eq. 57 to 
59 and rearrang 
the 
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_ illustrate the use of the modified three-moment equation in solving con- - 
tinuous” beam-columns on elastic foundations, the numerical example of Fig. 7 
-2(a) presented by the authors will now be considered, 
—— — of Eq. 60 to to spans ab and be leads to af 


a4 0, 29530 - 5.30923 + 2.25228 x 107° 


in which K = - 
_ To obtain the ‘second equation relating M,, Mp an and Yo» it is convenient to 


is assumed to be rigidly supported ai and infinite flexural rigidity. 
modified three-moment equation written for the imaginary —_ and the | span 


zt 0. 03671 — + 0, 01004 + 25228 x 10% 


62. 26324 x (64) 


ields: 


= 39294 4x 10 


‘These values arethe same asthe results obtained by the in 1 
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‘The equilibrium condition of Eq. 29(b) 
- 0.12034 + 4,60893 + 130,49145 
— 
Simultaneous solution of Eq. 62, 63, and 64 y | 
_M, = -19.5584 x 10 °K 
“My = 1.8688 x 10K 
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